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Abstract. In this paper we shall discuss the backreaction of a massive quantum scalar field on the 
curvature, the latter treated as a classical field. Furthermore, we shall deal with this problem in the 
realm of cosmological spacetimes by analyzing the Einstein equations in a semiclassical fashion. More 
precisely, we shall show that, at least on small intervals of time, solutions for this interacting system exist. 
This result will be achieved providing an iteration scheme and showing that the series, obtained starting 
from the massless solution, converges in the appropriate Banach space. The quantum states with good 
ultraviolet behavior (Hadamard property), used in order to obtain the backreaction, will be completely 
determined by their form on the initial surface if chosen to be lightlike. Furthermore, on small intervals 
of time, they do not influence the behavior of the exact solution. On large intervals of time the situation 
is more complicated but, if the spacetime is expanding, we shall show that the end-point of the evolution 
does not depend strongly on the quantum state, because, in this limit, the expectation values of the 
matter fields responsible for the backreaction do not depend on the particular homogeneous Hadamard 
state at all. Finally, we shall comment on the interpretation of the semiclassical Einstein equations for 
this kind of problems. Although the fluctuations of the expectation values of pointlike flelds diverge, if 
the spacetime and the quantum state have a large spatial symmetry and if we consider the smeared fields 
on regions of large spatial volume, they tend to vanish. Assuming this point of view the semiclassical 
Einstein equations become more reliable. 
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1 Introduction 

Despite the presence of many attempts, no widely accepted description of the quantization 
procedure necessary to treat quantum gravity is, up to now, available; however, nowadays we 
know how such a theory should look like, at least in some regimes. For example, treating fields 
as quantum object over a fixed curved classical background, many interesting effects have been 
discovered, like for example Hawking radiation in the case of black hole physics |Ha751 IFH90| 
or particle creation in the contest of cosmological spacetime |Pa68l IPa69j . Another interesting 
regime, where some progress could be made towards an eventual quantum gravity, is the analysis 
of the backreaction of matter quantum fields on curvature. In the present paper we would like 
to analyze some aspects of this problem in the context of cosmological spactimes. 

In this respect the recent developments in the theory of quantum fields on curved spacetimes 
provided by the algebraic approach to quantization |Ha92| . have been really fruitful. Here we 
think for example at the work [BFV03j where it has been shown, that the quantization is nothing 
but a functor that associates to every spacetime its corresponding algebra of fields/observables. 
It is thus possible to quantize simultaneously on every spacetime and, as a byproduct, the 
renormalization ambiguities, present in the definition of pointwise product of fields, acquire a 
clearer interpretation [HWOlj . In order to use similar ideas, a good control on the form of 
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the quantum states has to be available and, in this respect, an important contribution is the 
one of Radzikowksi |Ra96] (see also |BFK95j ). In that paper he has shown that, considering 
a free quantum theory and employing techniques proper of the Hormander microlocal analysis 
[Ho89j ■ it is sufficient to prescribe the form of the wave front set of the two-point function of 
a state in order to fix the form of its singularity (Hadamard form). Using similar ideas and 
similar methods it is also possible to show that many of the features of quantum fields on flat 
spacetime can be translated in a conceptual clear way to the curved case. More precisely, the 
Wick polynomial, their time ordered products and the stress tensor have been addressed in 
[HWOH IHW021 IHW051 IMo03] , whereas the study of the perturbative construction of interacting 
field theories can be found in [BFOnj \RWQ^\ iBFMj [BPFOQ] . 

Coming back to the analysis of the backreaction of quantum matter on classical curvature, 
a milestone is provided by the seminal paper of Wald [Wa77] . In that work, he has given a set 
of five axioms that a renormalization scheme should satisfy in order to produce a stress tensor 
whose expectation values can be safely used in the semiclassical Einstein equations 



in such a way that a physically reasonable backreaction emerges. Nowadays, the first four 
requirements on the expectation values of the stress tensor (consistency, covariance, conservation 
of {T^u) and standard flat spacetime limit) are automatically obtained in every locally covariant 
theory constructed as in |BFV03] (see also the nice lectures [BF09j ). However, the fifth one, 
which requires that the expectation values of the stress tensor depend on the derivatives of the 
coefficients of the metric up to the second order is more problematic. Nevertheless, even if it 
cannot be satisfied in general, it is possible to fix part of the renormalization freedom in order 
to have that at least the expectation values of the trace (T) := (T^^) meet this requiremeant 



Making use of these methods, in the case of cosmological spacetimes, the backreaction of a 
massless conformally coupled scalar field has been studied by Wald in |Wa78a| and in this case 
an exact solution can be found. Later on, another important contribution in this direction is the 
one due to Starobinsky |St80| (see also |Vi85j ) . where analyzing the trace of the stress tensor, and 
its anomalous term [Wa78b| . the backreaction on curvature was computed. Following similar 
ideas the analysis of the backreaction of massive scalar fields has been analyzed in |SS()2' and 
more recently in [DFP08] . Unfortunately while in the case of massless conformally coupled 
fields the backreaction can be analytically computed, essentially because the quantum state 
does not influence the dynamics of the spactime, when the fields have a mass, the situation 
is more complicated and the details of the quantum state have to be taken into account. A 
nice study of the semiclassical problem, also in terms of numerical integration, can be found 
in the series of works of Anderson |An85t IAn86j . The analysis of the problem of definition of 
the expectation values of the stress tensor in relation with the Casimir effect is presented in 
|GNW09] . The semiclassical Einstein equations have also been employed by Flanagan and Wald 
in [FW96] in order to study the validity of the averaged null energy condition on semiclassical 
solutions. Other interesting works at the interface of classical gravity and quantum matter 



are |GM86[ IRV991 lPS93l IPR991 lAEOOl IShOS] . Recently the dynamical system associated with 




(1) 



[Wa78b| . 
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the semiclassical backreaction has been studied in [EGlOj . where similar equations, as the one 
studied in the present paper, have been derived also for more general kinds of matter fields (with 
different curvature coupling). 

In this paper we shall follow the discussion presented in |DFP08| where some approximated 
solutions of the semiclassical Einstein equations have been displayed. In the latter work the 
matter has been described by a free scalar field whose coupling constant to the curvature is 
equal to 1/6 and with a non vanishing mass m. Notice that, despite the simplicity of the treated 
model, the solutions presented in that paper have very nice physical properties. On the one 
hand they show, at late times, a phase of exponential expansion (de Sitter), whose acceleration 
corresponds to a renormalization freedom of such a model, while on the other hand, at the 
beginning of the model a phase typical of the power law inflationary scenario emerges. In a 
certain sense that model provides an effective cosmological constant that varies in time and that 
could be used in order to model the present acceleration of our universe. 

Here, we shall show that exact solutions of that model exist and that many of the previously 
discussed physical properties shown by the approximated solutions persist. To this end, since we 
are here interested in the cosmological scenario we shall consider the categor}0 3" of spacetimes 
with a large spatial symmetry namely the flat Friedmann Robertson Walker (FRW). Thanks to 
the result stated in |BFV03j and [HWOlj . the fields realizing the components of the stress tensor 
can be unambiguously found on every spacetime, hence also on every element of 3". However, in 
order to have a well defined system of equations governing the backreaction, we must be able to 
unambiguously select a class of states (one for every element of 3") with good ultraviolet behavior 
in a spacetime independent way. Eventually, the validity of the semiclassical Einstein equations 
will be employed as a constraint on the class 3". The difficult part in this program is to give, for 
every element of 3", a state which is not based on the particular spacetime. Here we would like 
to give a prescription for obtaining those natural states. Furthermore, in order to give meaning 
to observables like the stress tensor, the states must have good ultraviolet behavior, namely the 
singular part of their two-point function must be as close as possible to the singularity shown 
by the vacuum in Minkowksi spacetime. Unfortunately this is not always possible and in many 
cases, a state, determined by some initial values on a spacelike Cauchy surface, does not satisfy 
the necessary regularity conditions. From the physical point of view, a natural candidate for a 
quantum state would be the choice of a vacuum, but unfortunately such concept is not available 
in curved spacetime. In any case we could ask for the validity of some physical properties in 
order to select a class of states that behaves like the vacuum in Minkowski spacetime, or at least 
that are close to it. To this end, Parker |Pa69| introduced the adiabatic states, namely states 
that are close to the vacuum. Afterwards, Liiders and Roberts [LJ190| made their definition 
precise employing the construction of Parker to select certain initial conditions at finite time 
in a FRW spacetime. Unfortunately, in general, the states obtained in that way do not satisfy 
the Hadamard condition, hence they do not have the same ultraviolet singularity as the vacuum 
in Minkowski spacetime, as shown by Junker and Schrohe [JS02]; nevertheless they satisfy a 

^The elements of 3" are the oriented and time oriented fiat FRW spacetimes and tlie morpiiisms can be cliosen 
to be the orientation preserving conformal transformations 
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condition weaker than the microlocal spectral condition given in |BFK95j . Another possibihty 
would be to consider states that minimize the smeared energy density. Olbermann |Q107| has 
constructed such states and has shown that they are of Hadamard form. Unfortunately, since 
this procedure depends upon the smearing over a finite interval of time, it is not possible to get 
those states only out of the initial conditions. 

In the present work we shall follow another procedure close to the one presented in |DMPn61 
IMo061 IM0O8] in the context of asymptotically flat spacetimes and already employed in cosmo- 
logical scenarios in |DMP09at IDMP09bj . We shall see that, if the initial surface is chosen to be 
lightlike (coinciding with the null past infinity in the conformally related spacetime) , we can give 
a class of states (that do not rely on the particular spacetime) that can be used in order to solve 
the coupled semiclassical system. Eventually, we shall comment on the late time behavior of the 
solutions showing that, if a de Sitter phase of exponentially accelerated expansion has emerged, 
the value of that acceleration does not depend on the particular homogeneous Hadamard state. 
Finally we shall give a new interpretation of the semiclassical Einstein equations when consid- 
ered as a non local equation. We shall actually show that, when smeared on large spatial volume 
regions, while the expectation value of the trace, (T), does not differ from the pointlike one, its 
fluctuations tend to vanish. 

In the subsequent section we shall introduce the geometric setup and the basic facts about 
the quantization we shall use in the next, this section is used to fix the notation and to recollect 
some results used throughout the subsequent part of the work. In the third section a discussion 
about the regularity and the Hadamard property shown by a class of states will be discussed, 
the states constructed in a spacetime independent way are presented in this section. The fourth 
section contains the main theorems about the existence of solutions of the semiclassical Einstein 
equations in a cosmological scenario. In the fifth section we shall discuss the behavior of the 
semiclassical solutions at late times and we shall see that the variance of the trace vanish when 
smeared on large volume regions. Finally, some comments and remarks are collected in the 
sixth section. The appendix contains a discussion about deformation quantization, the technical 
estimates and a technical proof of a proposition used in the main text. 

2 Geometry of the spacetime, classical and quantum matter 
2.1 Geometry 

Since we are dealing with problems in cosmology we shall assume homogeneity and isotropy. 
Hence, every element {M,g) in the class of spacetimes under investigation is such that the 
smooth manifold M is of the form M = I x T, where / C M is an interval representing the 
cosmological time and S is a three dimensional hypersurface that can be closed open or flat. 
The metric g, takes the well known Friedmann-Robertson- Walker form: 
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where dS'^{9,Lp) is the standard metric of the unit sphere in spherical coordinates and k G 
{+1,0, —1} distinguishes between the topologies of S, in particular between open, flat or closed 
spacelike hypersurfaces respectively. Furthermore, a{t) is the scaling factor representing the 
"story" of the universe and it is the single dynamical degree of freedom present in the system. 

In the next, also supported by recent observations that seem to confirm this hypothesis, we 
shall restrict our attention to the flat case, k = 0. With this assumption the metric ([2]) becomes 
conformally flat. This last fact, when a{t) is given, can be manifestly seen passing from the 
cosmological time t to the conformal time r. The latter is determined by the diffeomorphism 

where tq and to are two flxed constants. Since this transformation is a diffeomorphism and a is 
a field on the manifold M, it is clear that a can be seen both as a function of t or, as a function 
of r. For this reason, in the next, we shall indicate a o r~^(r) by a(r) as it is custom for scalar 
fields over manifolds. 

We would like to stress that choosing a{t) as the single degree of freedom, we are at the same 
time selecting a fixed (free) background on which a{t) evolves. The matter evolution could be 
analyzed on that background too, however, since we restrict our attention to conformally flat 
spacetimes, another natural choice for the background would be Minkowksi spacetime (M,(7m) 
taken with the standard Minkowksi metric 

gm ■■= -dr"^ + dr"^ + r'^dS'^iO, ip) , 

where (r, 6, ip) represents a point in in the standard spherical coordinates. In order to fix the 
notation, in the next, we shall indicate the generic element x of M by (t, x); notice that this 
defines a coordinate system also on M. Later it will also be useful, to consider as a background, 
the static Einstein universe {Me,ge). The advantage of this last selection is in the compactness 
of its spatial sections. All these spacetimes are connected by a chain of conformal embeddings 
[Pi09] rj and r]' such that 

{M,g)^{M,gM)-^iMe,ge) (3) 

where the conformal factor of rj is := a^^, while the one of rj' is 

^,2 ,^ |^ (l + (r + r)^)(l + (r-r)^) 

We would like to stress that in every spacetime under investigation, namely for every a{t), there 
is a conformal Killing vector, dr generating conformal translations r — >• r + A. We shall use 
these conformal symmetry in order to define a suitable quantum state used for the computation 
of backreaction. 
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2.2 Scalar fields, quantization and the stress tensor 

As a model for matter we shall consider a real scalar field, namely a field satisfying 



(4) 



For simplicity, we shall choose the coupling to the metric to be ^ = 1/6, corresponding to the 
so-called conformal coupling, whereas we shall allow the mass to be non vanishing. In this 
respect m shall play the role of the "coupling constant" for the matter gravity interaction, or 
the interaction between the field (p and the scaling factor a. 

We shall employ the quantization prescription proper of the algebraic approach ( |BFV03] 
IHWOlj ). First of all we select the *— algebra of local observables / fields (M) and afterwards, 
choosing a suitable state u (a normalized, positive linear functional over (M)), the expectation 
values for the elements of £/{M) are obtained. In order to keep the paper self-consistent, an 
explicit construction of (M) will be given in the appendix lA.il For our purpose it is important 
to stress that, being M a globally hyperbolic spacetime, {M) is uniquely determined by 
the equation of motion and the unique causal propagatocl E. The whole procedure is hence 
functorial as discussed in |BFV03] . In order to compute the backreaction of matter on gravity 
by means of the semiclassical Einstein equations ([T|), the expectation values of the component of 
the stress tensor need to be computed. Hence, since the pointwise product of fields is singular, 
the quantum state has to be selected in a physically reasonable way. In other words a certain 
regularization procedure has to be implemented and it must be meaningful on that state. 

We proceed now to discuss this aspect in detail. On the field algebra (M) (constructed as 
discussed on the appendix lA.ip . the functional u representing the state can be thought as being 
completely described by a suitable set of distributions a;„ € D'(M"), the so called n— point 
functions. In this paper we shall restrict our attention to states that are quasi-free, namely 
states whose n— point functions vanish for odd n while for even n they can be given in terms of 
the two-point function uo := L02 only by means of the formula 



where the sum is taken over the set P„ of permutations of the first n natural numbers with the 
constraints 

7r„(2z - 1) < 7r„(2i), 1 < i < n/2 and 7r„(2« - 1) < 7r„(2i + 1) , 1 < i < n/2 . 

Let us proceed to discuss the regularization prescription we shall use to give meaning to observ- 
ables like the component of T^y or the field square. In order to implement it in practice, we 
must have control on the ultraviolet behavior of our states. For this reason, we shall impose a 
constraint on the states, asking that their two-point functions have singularities that look, as 
much as possible, similar to the ultraviolet singularity displayed by the Minkowski vacuum. This 
is achieved requiring that the two-point function uj satisfies the so called microlocal spectrum 




(5) 



condition |Ra96[ IBFK95] . 



^The construction, and the proof of existence and uniqueness can be found for example in [BGP96| . 
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Definition 2.1. A state u: on s^{M) satisfies the microlocal spectrum condition (fiSC) if its 
two-point function to is a distribution on Cg°(M^) whose wave front set has the form 

WF{u:) = {{x,y,K,ky) eT*M^\{Q} , ix,k^) ^ {y,-ky) , k^>0} . 

The relation (x, k^) ~ {y, —ky) is satisfied in M if tlie points x and y are connected by a null 
geodesic whose cotangent vector in x is kx and in y is —ky. Furthermore /c^. t> if k^. is future 
directed. Notice that, thanks to the work of Sanders |Sa09j . it suffices to give a condition only 
on the two-point function. In the above mentioned paper it is actually shown that every state of 
a free scalar field algebra, whose two-point function satisfies the microlocal spectrum condition, 
has n— point functions that enjoy the microlocal spectrum condition too (for the definition in 
that case we refer to the work |BFK95| ). 

Notice that, the two-point function of a quasi- free state for fields satisfying the equation of 
motion Q and the microlocal spectrum condition, has an ultraviolet singularity that is similar 
to the one of the Minkowski vacuum. In other words, its integral kernel can be expanded in the 
following way 

a;(x, y) := lim f^lpA + v{x, y) log + W{x, y) (6) 

where U, V and W are smooth functions on , A is a length scale used in order to have the 
logarithm with a dimensionless argument and a^{x,y) = a{x,y) + ie{^{x) — -^{y)) + where 
^ is a generic time function and a is the Synge world function, namely half of the squared 
geodesic distance taken with sign [D BGOj iKWQlj . Furthermore, e is a regularization parameter 
that eventually tends to 0"*" (weak limit). A two-point function whose singular part is of the 
form ([6]) is said to be of Hadamard form. It is important to notice that the coefficient U 
and V are completely determined out of the geometry and of the equation of motion, while the 
coefficient W really characterizes the state. The distribution constructed out of the first two 
terms in ([6]) is called Hadamard singularity and it is indicated by Ji, then its integral kernel 
is 

JC(x, y) := hm f + V{x, y) log . (7) 

Hence, in order to regularize the states satisfying the microlocal spectrum condition we can 
simply subtract the singular part "K from cj, what remains is a smooth function whose coinciding 
point limit can be safely computed. As a side remark, we notice that, since we are interested in 
local fields, for our purpose it is enough to have an explicit form of "K on small domains. The 
question about global existence of "K is hence not an important one in our approach. Coming 
back to the main point, it is hence possible to obtain a well posed stress tensor on states satisfying 
the /xSC. On practical ground this is achieved by operating on the regularized two-point function 
Lo — 'K with a certain differential operator before computing the well posed coinciding point limit. 
This procedure is called point splitting regularization. Furthermore, we notice that on those 
states the fluctuations of the observables computed on similar states are always finite. 
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Let us discuss this point in some details. On the class of states that satisfy the fiSC, we could 
extend the algebra of fields (M) in order to encompass also pointwise product of fields. Here 
we shall indicate by ^(M), the extended algebra of fields, namely the *— algebra that contains 
also the Wick powers as constructed in [HWOlJ . or more recently by means of deformation 
quantization methods in |BDF091 IBFOQj . For completeness, the main steps of that construction 
can be found in the appendix EH We are now ready to discuss the form of the stress tensor as 
an element of ^(M). To this end, following the discussion introduced in |Mo03j . we consider 
the following operation on smooth field configurations (p 

Tab{(p) ■■= dafdbif - ^gab {dcfd^lf + m^(/?^) - adb'/ + i (^Rab - ^9ab^ V?^ + - gab^'p'^- 

Notice that for every compactly supported smooth function /, 

Tabim) ■.= l{f,Tab{v>)) 

is an element of ^(M), and we shall use it in order to compute the backreaction in ([T|). Another 
element of ^{M) that will be useful is 4>^{f)- This is nothing but J ^ifp' jl. We shall indicate 
by (<^^)a; and {Tab)uj the expectation values of the two observables introduced above computed 
on the state u regularized according to the previous discussion. Hence {4>^)u: and {Tab)uj are 
obtained replacing the smooth field configuration ip{x)(p{y) with the smooth function oo — "K 
in Tab{^){f) and (jP'if) before performing the coinciding point limit (contracting with the delta 
function). Further details on this regularization, can be found in appendix lA.il In the procedure 
presented above, there is an ambiguity in the definition of IK, whose form is described in |HW01j : 
we call it renormalization freedom, later on we shall give more details on this crucial aspect. 

3 Homogeneous states and Hadamard property in flat FRW 

3.1 Characterization of pure, homogeneous states in terms of solutions of 
the wave equation. 

In this subsection we shall proceed to discuss the form of the states we are interested in, we 
shall make use of some previous works in order to get the form of the two-point function of 
the pure homogeneous states on FRW spacetimes in particular. We start with the definition 
of the class ^{M) of pure and homogeneous quasi-free states, we shall employ in the present 
paper. Although the set "^(M) does not contain all the possible pure and homogeneous states 
for =2/(M), that definition is motivated by the works of Liiders and Roberts |LR90] (Theorem 
2.3 of | LR90) in particular) and the one of Olbermann |O107j see also the nice work of Degner 
and Verch |DV09j for a review. Furthermore, we shall restrict our attention to that class. 

Definition 3.1. We shall indicate by '^{M) the set of pure, homogeneous and quasi-free 
states u) for j/(M) whose two-point function to S D'[M x M), using the standard Minkowksian 
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coordinates, takes the following form, in the distributional sens^, 

a;((r,x); (r',x')) := 7^ / ^^e^'^-(--)<ik , (8) 

where (r, x) is a generic element of {I xW^ , g) , k = \k.\. For every k > 0, Sk is a smooth function 
satisfying 

^Skir) + {m'airf + k^)Sk{T) = 0, (9) 



and 



Sk{r)^Sk{T) - ^Sk{r)Sk{r) = i . (10) 
(It (It 



Furthermore, for every r, the map k i— t- Sk{T), and its first time derivative, are measurable 
functions that are polynomially bounded at large k. When restricted on a generic interval [0, /cq] C 
M+ they are in L^{[0,ko],k^dk). 

It is not guaranteed that the state, constructed out of a generic solution Sk of the equation 
([9]) (also with the constraint on the Wronskian), is of Hadamard type. It is also usually not 
easy to impose some ad hoc initial conditions for S'fc(r) on a spacelike hypersurface in order 
to have an Hadamard state. This problem, in connection with the concept of adiabatic states, 
has been discussed in the paper of Junker and Schrohe |JS02j who showed that it is possible to 
get Hadamard states only if the recursive construction of adiabatic vacua does not break down. 
However, as shown by Liiders and Roberts |LR90] . one cannot exclude such possibility. Another 
possible explicit construction of Hadamard states is the one presented by Olbermann |O107j . 
but this construction requires the smearing on a finite time interval. It is hence very difficult to 
use similar ideas in treating backreaction problems. Here we would like to give an alternative 
construction employing ideas similar to the one presented in |Mo08l IDMPOOal IDMPOObj and 
adjusting them in order to be suitable for the purpose of studying the gravity-matter interacting 
system. First of all, we start selecting a particular solution Xki^) of ([9]), when the spacetime M 
admits the limit r — )• — oo, namely the one that satisfies the initial conditions 

hm e'>'-xk{r) = ^, lim e^^^ ^Xkir) = . (11) 

r-^-oo ^2k r-*~oo (iT V Z 

We can now introduce the following proposition that summarizes some results obtained in [LR901 
[QlOTl IDV09I 

Proposition 3.1. If the limits in ^21\) are well posed, every solution Sk of (0) which determines 
a state u: G ^(M) according to the definition \3.1\ can be realized as a linear combination of the 
Xk (the solution of the equation (0) satisfying the asymptotic conditions (fii])).- 

Sfe(T) := A{k)xk{T) + B{k)xk{T), (12) 



'Here we mean that ui has to be constructed as the weak hmit e — )■ 0+ of the expression on the right hand side 



of ((8]) regularized multiplying the integrand by an extra factor 



e 



-elk 
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and the coefficient A(k) and B(k) are such that 



\Aik)\'-\B{k)\' = l. 

Furthermore, for every t, the functions, on R^, k i— )• A{k)xk{'T) and k i— )• B{k)xk{'^), together 
with k I— )• A{k)^Xk{'T) and k i— )• B{k)-^Xk{T), are at most of polynomial growth and locally 
square integrahle. 

Notice that, even if we restrict our attention to the class of pure stated, we have still some 
freedom in the definition of the state, as expressed in the preceding proposition by the choice of 
A{k) and B{k) in the Sk of ()12p used in finding the two-point function ([8]). Motivated by this 
expression, we shall indicate by u>a.b the quasi-free state whose two-point function is 

T,x ; r ,x ) := — --3 / dk e ^ > —- -— (13) 

(2vr)'^ J^i aij) a{T') 

where we have highlighted the dependance upon A{k) and B{k) in the state. 

We are now going to find sufficient conditions that have to be fulfilled by A{k) and B{k) 
in order to ensure that uja.b G ^(M) satisfies the Hadamard condition. The strategy we shall 
pursue is the following one. We prove that ujifi is of Hadamard form provided some regularity 
is shown by a{t) in ([2]) and then that u:a,b satisfies the /xSC if and only if B(k) is rapidly 
decreasing for large k. Notice that, in order to prove that a certain distribution satisfies the 
microlocal spectrum condition, it is not enough to test it on some Cauchy surface because no 
remnant of causality survives the projection of any vectors on a spacelike hypersurface. This 
is not true for projection on lightlike hypersurfaces, and such is the reason why the arguments 
used in [MoOSj . [Honn] and |DMPn9bj works. We would hke to give sufficient conditions on 
the form of the spacetime, and on its past asymptotic form in particular, in order to guarantee 
that ui^o is of Hadamard form. In the case of asymptotic de Sitter spacetime, it was already 
shown in [DMP09bj that the state u^ifi satisfies the microlocal spectrum condition, and that 
if the spacetime is precisely the de Sitter one, oji^q becomes the well known Bunch Davies 
state |BD78| . Furthermore, in the case of asymptotic flat spacetimes, Moretti has proven the 
Hadamard property for an analogue state [M0O8] : a similar proof can also be found in [HoOOj . We 
are now going to generalize these results in order to encompass the case considered here. The idea 
we would like to use, is to consider uJifi as the pullback (rjor]')* (uj) associated with the conformal 
embedding rj o rj' introduced in ([3]) where u> is the two-point function of a certain state in a 
hyperbolic subset of the Einstein universe that contains the image of M under r/'ory together with 
its past boundary Ui~. Notice that on the latter, the field theory looks like a massless Klein 
Gordon field perturbed with an external potential of the form m^a(r)^(l+(r+r)^)(l+(r— r)^)/4, 
and hence a) is a state for this particular theory. It happens that, if that potential becomes 
regular on the extended spacetime in a neighborhood of $5" U i~ then the state <^i,0) defined in 
such a way, turns out to be of Hadamard form. 

^More details in this respect can be found in |LR90l [O107) 
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Theorem 3.1. Suppose that M = / x R'^ where I = (0,to) (^f^d that a{t) is such that the region 
t — )• corresponds to the past null infinity of (M, (^m) in (^ejS'e) md furthermore such that the 
smooth function 



(l + (r + r)2)(l + (r-r)2)m2a(r)2 



(14) 



defined on M seen as embedded in the Einstein universe {Me,ge) can be smoothly extended in a 
neighborhood C Mg containing the past boundary of M here indicated as Q^Ui^. Under these 
hypotheses the state u^ifi, constructed out of the solutions Xk that satisfy the initial condition on 
Q~ satisfies the fiSC. 

Proof. The proof can be found following similar ideas as those presented in |Mo081 IHoOO) or 
alternatively in jDMPOObj . we shall summarize here the main strategy and we shall discuss only 
the key point in detail. We would like to show that the two-point function loi^q satisfies the 
microlocal spectrum condition, hence we have to analyze its wave front set and to show that 



Notice that, in showing the preceding equality, the difficult part is to prove the inclusion C. If 
it holds true, the other inclusion descends from the fact that uifi is a state |SVOH [SVW02j and 
out of an application of the Hormander's theorem of propagation of singularity jDH72j . In order 
to show the validity of inclusion C we can rewrite 0Jifi{f,g) as Ki{QEf \c^-,QEg |"cj_) where to 
is a certain distribution on 9^ x 3^ defined employing the Fourier Plancherel transfornj^ along 
the T— coordinate on 



J§2 Jo 

where is the measure on the two dimensional unit sphere and ^ is a shortcut for the 
standard spherical coordinates. A discussion on such a distribution can be found in theorem 
3.1 of |Mo08j . Furthermore, i^Ef \ 9~ is the restriction to 9~ of the wave function E{f) 
constructed by means of the causal propagator E, multiplied by a certain factor 0,. That 
restriction can be safely computed in an Einstein universe and the multiplication by the factor 
O corresponds to the conformal transformations that maps solutions of P in M to solutions of 
in Mg. We have indicated by Pe the operator -De + ^Re + rn^a^^'^ where is the d'Alembert 
operator associated with the Einstein metric. Hence, the two-point function cji_o can be seen as 
a composition of distributions 



Because of the non compact nature of the restriction of E on 9", which persists also if the other 
entry is localized on a compact set of M, the previous composition of distributions cannot be 
straightforwardly justified by an application of theorem 8.2.13 of |Ho89) . On the other hand, in 

^This operation is introduced and discussed in the appendix of [Mo08] . 



WF{uJifi) = {{x,y,k^,ky) £T*M^\{Q] , {x,k^) ^ {y,-ky) , A;^>0} 



(15) 




uifi :=tvo {QE fcj- ^QE fcj-) . 
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proposition 4.3 in |Mo08] . it is proven that a similar composition is well defined and that the 
result has the desired wave front set provided the following two facts hold: First of all considering 
a compact set inside M, there exists a neighborhood Oj- of i~ in Me such that there are no 
lightlike geodesies connecting any point in with those in Oj-, and, second, if one entry of 
the causal propagator is restricted on and the other on 9", the causal propagator itself must 
fall off sufficiently fast towards i~ in M. Both requirements hold under the hypotheses of the 
present theorem, the first descends from the fact that M is conformally flat and the latter is a 
consequence of the fact that m'^a{Tf{l + {T + r)^){l + {T-r)^) can be smoothly extended on a 
neighborhood of in Mg, hence the same result as the one stated in lemma 4.2 of |Mo08| holds 
in the present case too. □ 

Notice that the constraints given above are only sufficient but not necessary, in fact in |DMP09bj 
it has been proven that cji^o is a well defined Hadamard state also when the spacetime is 
asymptotically de Sitter in the past, and in such a case the potential that appear in the equation 
of motion on the Einstein universe is singular on i~ . Hence the hypothesis of smooth extension 
of the function (1140 over 3" U i~ in Mg can be relaxed. Finally, notice that employing similar 
techniques to the one presented in [DMP09a] and in |DMP09b] based on a careful analysis of 
the modes Xk{T), it is actually possible to obtain the Hadamard property for the state uifi 
also when (|14p can be extended over 9~ U i~ only with the regularity of C^(Me). We would 
like to conclude this section by introducing a theorem that gives some necessary and sufficient 
condition for the state uja,b to be of Hadamard form. 

Theorem 3.2. Suppose that coi^, constructed as in (jl3p . is of Hadamard form. Then the 
quasi-free state ua.b £ '^{^), satisfies the microlocal spectrum condition if and only if |-B(A;)| 
decreases rapidly for large k, in other words, there are constants c„ and k such that 

\k\''\B{k)\ <cn , y\k\>k. 

Proof. The proof of this fact can be easily done using the formulation of the state on a given 
Cauchy surface, see |LR90j for a discussion. Notice that, when restricted on a Cauchy surface, 
the difference between the two-point functions u}a,b — ^i,o becomes smooth if and only if B{k) 
decreases rapidly (and this holds for the time derivative of one or both of its entries too). This 
proves that the rapid decrease of |-B(fc)| is a necessary condition for the validity of /xSC. Finally, 
in order to prove that it is also sufficient, we can make use of a formula like the one given in 
equation 3.60 of |KW91j and reconstruct uja,b — ^i,o out of its form on a Cauchy surface. We can 
extend the result about the smoothness to the whole by noticing that every of the four terms 
appearing in the formula 3.60, above cited, can be obtained composing u}a,b — ^i,o (sometime 
with a time derivative applied to one or both of its entries) with the tensor product E ® E 
(sometime taken with a time derivative applied to one or both of its entries). The proof can 
be easily finished noticing that E maps compactly supported smooth functions on the Cauchy 
surface to smooth functions in the spacetime. □ 
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4 Semiclassical treatment of backreaction and existence of so- 
lutions 



4.1 Semiclassical approximation 

We shall now proceed to treat the backreaction of quantum matter on gravity, and, as discussed 
in the introduction, we shall do it for flat FRW spacetimes. With such a large symmetry, the 
Einstein equations take the well known Friedmann Robertson Walker form, namely, being the 
stress tensor of matter diagonal and equal to T^'^ = diag{—p,P,P,P), if we indicate by a the 
derivative of a with respect to the cosmological time t and H := a~^d they become 

3/7 = Svrp - -2 , 3(^H + H^j = -47r(p + 3P) . 

Because of the conservation of the stress tensor we can instead solve the equation for the trace 
of Gf,u and of T^u 

-R = SttT 

and then employ the first FRW equation as a constraint, that needs to be imposed only at a 
fixed time, in order to fix only an initial condition. We shall assume this point of view and, in the 
semiclassical scenario, we have simply to substitute the classical T with the expectation value 
of the corresponding observable (T) computed in the chosen quantum state. More precisely, 
as discussed in the introduction, we have to select a class of states (one for every a{t)) and 
then the semiclassical equation —R = 87r(T) acts as a constraint on this class. Notice that, by 
employing the initial values on the null boundary, we have already disentangled the problem of 
the definition of the class of states from the geometric property of the spacetime. 

Hence the backreaction of such system can be treated simply employing the equation for the 
trace. But in the trace, taking the conformal coupling and a non-vanishing mass, the explicit 
form of the state become important. More precisely, the expectation value of the trace of the 
stress tensor in a state u: is 

Above [Vi] is the coinciding point limit of the second Hadamard coefficient (see [SVOlt IMo031 
IDFP08| for more details). This contribution to the trace is the anomalous one, it coincides 
with the well known trace anomaly displayed by the stress tensor |Wa78bj in the massless case. 
Coming back to the dynamical equations, in the case of a fiat FRW metric, the equation of 
motion for the scaling factor is 

- 6 (i? + 2H^^ = -87rm2(</.2)^ - -1- (^HH^ + H^^ . (16) 

In writing the preceding equation we have performed some choices of the renormalization pa- 
rameters, in particular we have chosen some of them in such a way that the higher derivatives 
of the coefficients of the metric disappear. There is an extra renormalization freedom that we 
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have not considered here, but it can always be reabsorbed in the freedom present in It 
seems important to notice that, apart from these renormahzation freedom, the trace anomaly 
is fixed and it is really a c— number. For the discussion of these and other points we recall the 
work [DFPnSj . 

Hence, in the case of vanishing m, the state does not enter anymore in ()16p and an exact 
solution can be easily founci^. This solution can be written in an implicit form as 



H + H. 



H-H. 



(17) 



where the time of the big bang is fixed to be t = and = 3607r in natural units, that is 
G = h = c=l. 

The massless solution (I17p presents three branches for H > (the first < < the 
second < < H\ and the last > H'^), we concentrate ourself on the upper one and 

we discuss some of its properties. Notice that, it corresponds to a universe that shows, at the 
beginning, a phase typical of the power law inflationary scenario, while it ends up in a phase 
which corresponds to a de Sitter universe. It is a remarkable fact that the form of the initial 
singularity is lightlike; hence, for example, the horizon problem is not present in this simple 
model. The Hubble constant H of the de Sitter phase shown in the asymptotic future is equal 
to -fr+, and it is many order of magnitude too big in order to explain the present acceleration 
of the universe. Despite of this fact, in [DFPOSj, it has been shown that, in the case of massive 
fields, this parameter is not fixed and it needs to be interpreted as an extra renormahzation 
freedom. In the next we shall make use of the lightlike nature of the initial singularity to see if 
it is possible to obtain an exact solution also in the case of massive fields, and we shall look for 
solutions that are similar to the massless one at the beginning of the universe. In this respect, 
we shall treat as a renormahzation constant. Furthermore, since, in the massless case, only 
the anomaly appears in the trace and it is a c— number it has a vanishing variance, hence, such 
equation should hold also for large H (as is the case close to the initial singularity t = in the 
(jl7p ). Of course, the situation changes when the mass is different from zero, and in this case the 
result provided by the semiclassical Einstein equations needs to be interpreted more carefully. 
We shall comment on this points in the next section. 



4.2 Initial conditions at the beginning of the universe 

In this subsection we would like to discuss the existence of solutions of the semiclassical Einstein 
equations ([T6|) that look as much as possible close to (fT7|) . at least on small intervals of proper 
time just after the initial singularity. Hence, we shall assume the form of the singularity to be 
lightlike just as in the massless case (fT7|) . As previously discussed this task can be accomplished 
only after a careful analysis of the class of states we have to use in order to obtain the expectation 
values of (j)'^ in the different spacetimes we are considering, in any case the requirement of having 
a lightlike initial singularity is essential in order to apply the results of the preceding sections 
and to have Hadamard states. 

•'see |Wa78al IDFP08| for more details 
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The states we would like to chose are the one constructed in the preceding part of this paper 
and denoted by <^a,b with B of rapid decrease. We shall start considering only the state u^ifi, 
and we notice that these states are defined by the requirement to be vacuum states with respect 
to the conformal time on the initial singularity (this procedure works only for lightlike initial 
singularities) . Before proceeding with the discussion of equation (|16p we have to discuss the form 
of the expectation value of 0^ on these states and to fix once and forever the renormalization 
freedom, according to the picture presented in f BFV03| . 

Since we would like to impose initial conditions on 9^ or better at r — )■ —oo, where the 
Hubble constant H usually diverges, it is useful the rewrite the dynamical equation for its 
inverse that we shall call X := . With this in mind, we can rewrite the equation (116p as 
follows 

dt ~ 1 - i72 X2 + 1 - X2 '"^^'^ ' ^^^> 

where := H\/2 = ISOvr and := 2407r2m2 are two constants. Notice that, since {(j)'^)uiifi 
is a functional of a{t) and of r, the proof of the existence of solutions of (fT6|) needs some effort, 
and the analysis of the functional dependance of (0^)a)i_o upon X = needs to be carefully 
performed in particular. 

We intend to prove the existence of solutions of (116p (or better of (llSh ) in a certain compact 
set Sc of a Banach space To this end let us start introducing these spaces of functions. 

Definition 4.1. Consider the set ^ as the subset of elements of C^{[0,tQ]) that vanish in and 
whose first derivative is bounded on [0, to] . is a Banach space when equipped with the norm 

Wfl^ := sup fit) 
[0,io] 

For every | < c < 1 we define !Bc as the closed ball in of radius 1/2 — c/2 centered in 

/o(t):=i^t (tG[0,to]). (19) 

The strategy of the proof of existence of solutions of (jlSp is the following: We shall build a 
sequence X„ in !Bc by a recursive use of a certain map T : — >• defined as 



T(X) it) 



1 — X'^ Cm X 



1-mx^ i-mx^ 



dt' (20) 



and we shall show that, since 7 it is contractive on "Be, Xn converges to a solution. Hence the 
existence and uniqueness descends from an easy application of the Banach fixed point theorem. 
Notice that, if we succeed in obtaining that result, since certain physical properties (like the 
presence of a phase of power law inflation) are enjoyed by every element of 23c those properties 
have to be shown by the solution too. 

We shall now discuss how it is possible to associate a spacetime to an X G in such a way 
that the initial conditions are implemented. 
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Proposition 4.1. Fix the three constant to, ag and tq. Then every element X ofB^ determines 
a spacetime {M,g[X]) where M = [0, to] ^ o'^'^ where g[X] is the metric ([2]) constructed out 
of the following scaling factor 




■to 



) 



a[X]{t) := oq exp 



x{t')-^dt' 



The Hubble parameter shown by {M,g[X]) is then 



X{t) 



-1 



and the corresponding conformal time 



r[X]{t) :=ro 




a[X]{t')-^dt' 



is a functional of X too. Finally, the spacetime enjoys the following initial conditions 



X(0) = 



a{to) = ao 



T(to) = To 



(21) 



per construction. 

For every X in Sc we have immediately, from the definition, that the initial conditions ()21|) 
given above are satisfied, and the functional a[H] and t[H] are, moreover, well defined on 23c. 
Notice that, on Sc 



out of this, we have that a(0) = sufficiently fast in order to ensure that t[X] — t- — oo for t — t- 0. 
This last requirement is necessary in order for the resulting spacetime to have a lightlike past 
boundary or, in other words, that the initial surface t = is of null type. As already discussed, 
this last requirement is a necessary prerequisite for the construction of Hadamard states on 
(M,g). We have in fact, the following two propositions 

Proposition 4.2. For every spacetime characterized by a metric of FRW type, whose scaling 
factor a{t) is such that X{t) = is in "3^0 0°° {[0, to]), we can construct ojifi as in (jl3p and it 
determines a quasi-free state uji^q which, moreover, satisfies the microlocal spectrum condition. 

Proof. If the analogue of function (jl4p can be smoothly extended over S^Ui^ the proof descends 
from theorem 13.11 otherwise, also according to the discussion given after that theorem we can 
adapt the proof of theorem 3.1 given in pMP09b] to the present case. □ 

Notice that, if we drop the requirement of being smooth, and we replace it only with X £ 
"Be, the construction of the state is still valid, and, even if the resulting two-point function is 
not of Hadamard form, we can still apply the procedure to compute expectation values of (f^. 
Unfortunately, for fields that contain derivatives, it is not anymore guaranteed that the point 



lim 




= OO 
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splitting regularization works for these states (they are adiabatic states of suitable order, in the 
sense given in |JS02| ). 

We shall conclude this subsection by noticing that, since a[X] is positive, r[X](t) can always 
be inverted, furthermore, on Sc some useful inequalities hold. We shall summarize them in the 
next proposition: 

Proposition 4.3. The following inequalities hold for every X in 'B^. with 2/3 < c < 1 and for 

te[0,to] 

X(t) t 1 /t\- a\X](t) t 

t X(t) c V*o/ «o ^0 

The proof descends from a straightforward application of the definition of a and r out of H 
and from the bounds satisfied by the elements of Sc. 

4.3 Expectation value of (f? on the states oji o, and further renormahzation 
freedom. 

We are now in place to discuss the form {(t)'^)uii^o constructed out of point splitting regularization 
on the quasi- free state a;i^o introduced above in (fT3|) . namely we can construct it by 

{^^{x))u:^ = lim [^^1,0(3;, y) - ^{x, y)] + aR{x) + ^m^ (22) 

where IK is the Hadamard parametrix ([7]) constructed in the spacetime {M,g), and a and f3 
are renormahzation constants. At this point it is important to stress that a and /3 are really 
constants, in the sense that they do not depend on the form of the spacetime. In other words, if 
they can be fixed in a specific spacetime, they are fixed on all of them. This is a consequence of 
the general covariance presented in [BFV03j which uses the hypothesis that the fields transform 
covariantly under isometries used in combination with spacetime deformation arguments. 

According to the discussion presented in [ DFP08] we would like to have Minkowksi spacetime 
as a solution of the semiclassical equation (|16p . when the Minkowksi vacuum is employed. This is 
tantamount to require that the expectation value of 0^ on the Minkowksi vacuum in Minkowksi 
spacetime vanishes and this fixes the renormahzation constant /3 to a certain value. We notice 
also that the other renormahzation constant a can be used in order to adjust the parameter 
He present in the equation (jlSp . hence from now on we shall consider it as a renormahzation 
constant |DFP08j . 

Once we have a prescription to fix the renormahzation freedom we can start to discuss 
the Hadamard regularization; to this end we need to know the Hadamard parametrix '}{ in 
some details. For our purposes, since the modes x ^-re easily treated when analyzed on the 
conformal Minkowksi background, it would be desirable to perform the subtraction loi^ — IK 
directly on Minkowksi spacetime. To this avail we have to analyze how oji^q — '}i transforms 
under the push forward r/^, associated to the conformal transformation r] : {M,g) — )• (M,5fM) 
introduced in ([3]). Notice that r/^, maps uji^ to a(ti)a(t2)wi,o seen as a distribution on Minkowski 
spacetime, and, moreover, since the conformal transformation preserves the Hadamard property 
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|Pi09j ■ also a{ti)a{t2)u}ifl satisfies the /iSC in M. We can hence regularize it by means of the 
Hadamard parametrix "Km constructed in Minkowksi spacetime for the field theory Pm'/' = 
with Pm '■= — Dm + a^m?. Notice that, under the conformal transformation rj, the mass term 
on M becomes a position dependent potential a{t)'^m? . 

In this respect, we already know that the microlocal spectrum condition, together with the 
causal propagator, are preserved by a conformal transformation [Pi09] . and this holds also if the 
conformally coupled quantum field has a mass. Of course, under a conformal transformation, 
the symmetric part of the Hadamard function is not preserved, hence the coinciding point limit 
of the difference does not vanish in general, actually, taking the same scale A in the logarithmic 
divergence, we obtain that 



lim 



'}i{xi,X2) 



a{Ti)a{T2) 



'Km{xi,X2] 



m 



loga(r2) + cR{x2) 



where c is a fixed constant. Notice that an unexpected logarithmic dependance in the scaling 
factor log(a(T)) appears and we have to carefully use it. On the other end, the term proportional 
to R, which appears when the two different regularization schemes are used, falls into the class 
of the standard regularization freedom [HWOlj and, thanks to the previous discussion we can 
incorporate it in the definition of He and forget it. 

Furthermore, since we are only interested in computing the expectation value of c/)'^, it is 
possible to regularize by subtracting an approximated version of the Hadamard parametrix 
constructed just out of the first Hadamard coefficient (see [SVOH IMo03| for more details) 



rw-U 



lim 



According to the analysis performed in chapter 4 of |Fr75j . the Vo{x,x') for spacelike separated 
X and x' can be found as the result of the following integral 



Vo{x,x') := -Uix,x') 



f 


'PmU' 


Jo 


u 



{x,-f{s))ds 



where the integration is done in the affine parameter s of the geodesic 7 connecting x and 
x' with 7(0) = X and 7(1) = x' . The previous integral can be easily performed on Minkwoksi 
spacetime where U = l/(87r^), Pm = — DM+m-^a^, and the geodesies are simply the straight lines 
connecting x and x'. Furthermore, if we choose x and x' as (r, x) and (r, x') in the Minkowksian 
coordinates, we obtain 

Vo{x,x') = j^-^rn^aiTf . 

We can now write the Hadamard parametrix at first order, for x and x' lying on the same 
Cauchy surface S,-, as 



x') :- 
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Notice that the first two contributions are the same as the one of a massive Klein Gordon field 
propagating in Minkowski spacetime with the mass equal to m^a^(r) (where here r is fixed by 
the choice of E,-). Furthermore, we know that in this case, the coinciding point limit of the 
following expression yields 



lim 



'K^{x, x') 



1 



1 



(2^)3; 2^k2 + m2a(r)2 



gik(x-x') ^J^ 



rn?o? 



(4^) 



which has to be understood in the sense of distributions, while c is a constant. Hence, as a side 
remark, it is clear that the Hadamard regularization and the first order adiabatic regularization 
|Pa69] for the computation of the expectation values of cfP' coincides up to the choice of a 
renormalization constant. Furthermore, with a judicious choice of A (proportional to m), a 
direct computation shows that for every x and x' lying on the same Cauchy surface S,- 



lim 



"^mix, X 



(27r)2 



e(|k| -ma(r)) 



1 



2|k| 



m?a{Ty 
4|k|3 



gik(x-x') ^J^ 



(- log a+c') , 



and hence we could use the right hand side to regularize the states in order to obtain the 
expectation value of (fP' . Once again we incorporate the constant c' in the renormalization 
constant /3, that we shall fix a posteriori according to the discussion given above. 

Eventually, also in the computation of the expectation value of (f)^^ we are free to restrict 
both X and x' to the same Cauchy surface, determined by ri = T2 = t and, after the subtraction, 
we can perform the coinciding point limit 



(27r)3a(r)2 



(02(x'))^,,„ := 
Xfc(r)Xfc(r)-e(|k| 



ma{T)) 



2 k 



4|k|3 



(23) 



where /3 can now be fixed to the value — l/(87r^) by the requirement that cfP' vanishes in the limit 
of flat spacetime, furthermore, the integrals present above are to be intended in the distributional 
sense. Notice that the log(a) terms are not manifestly present anymore in the expectation value 
of (/>2. We conclude this section with a remark. Since the map t — )• r is a diffeomorphism, 
the previous expression can be easily given in cosmological coordinates, just employing a{t) 
and Xk{t) at the place of a(r) and Xki^)-: where, as discussed above, a(r) = (a o r~^)(r) and 
Xkit) := {xk°T){t). 



4.4 Recursive constructions and bounds satisfied by Xk and by the auxiliary 
functional $ 

The aim of the present section is to find suitable bounds satisfied by ui^Q^cp'^), for small values 
of a and close to t = 0. In order to accomplish this task we need some better control on the 
explicit form of the Xk appearing in ()23p . Hence, here we shall give a recursive construction 
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for Xk^ used in the definition of uji^q using m as a perturbation parameter. While the analysis 
of the state, and of the modes Xk in particular, is easier if performed in the conformal time r, 
the form of the differential equation (jl6p suggests the use of the cosmological time t. Here, as 
previously said since the map r from the cosmological time to the conformal is a diffeomorphism, 
without ambiguity we shall write Xk{i) for {xk ° where Xk{T) is a solution of dH) with the 
initial conditions (llip . Following similar ideas to the one presented in |FH90j though in another 
context, we can write XkiT) by means of perturbation theory, hence in term of a Dyson series 
over the massless solution Xk treating m?a{T)'^ as the perturbation potential (a similar 
construction can be found in the proof of Theorem 4.5 in |DMP09a] ). Finally, since we would 
like to use the cosmological time instead of the conformal one, we can simply compose Xk with 
r (writing Xki'T'i't))) and then change the coordinates in every integral appearing in that series 
from the conformal time to the cosmological one. Following this procedure we obtain 



Xfc[a,T](0 = X^X^[a,T](t), (24) 



n=0 



where. 



and 



Xlla.rm = -f^ MHr{t)^ ^^''^^K itWxr'[<^,r]it') dt' , (25) 
Xk[a,r]{t) = 



'2k 

The advantage of using the cosmological time is in the fact that it becomes clear how Xk depends 
on r. Notice that r appears only in sin(A;(r(t) — T{t'))) of ([25]) and in the form of x^- This 
decrease the number of functional derivatives that needs to be performed in the computation of 
the functional derivatives of Xk with respect to the small variations 5X of X. The first thing 
we would like to check is that the perturbative construction is well defined, namely let us show 
that the sum in (|24p converges uniformly. We shall state the result about convergence in the 
following proposition and, although the proof is almost straightforward we would like to give it 
since it contains some useful remarks that we shall use also later. 

Proposition 4.4. Let {M,g) be the flat FRW spacetime determined by X £ 'Be according to 
the statement of proposition \4.1\ where a is its scaling factor and r the corresponding conformal 
time. Then the series (f24]) . constructed over {M,g), converges absolutely to Xk on [0, to]; where 
Xk is a solution of equation ([9]), when read in conformal time, and, moreover, it enjoys the 
initial condition (jlip . Furthermore Xk' given in (j24p . and Xk satisfy the following bounds 

, r,/ M 1 (mt)'^^ , 1 fm?a{t)t\ , 1 / 9 9\ 



Introducing the following notation 

00 

Xn^,rm= Xl[a,r]{t) 

=n+l 
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we have also that 

Proof. By direct inspection, under the hypotheses stated in the proposition, if the series con- 
verges, it has to tend to a solution of the equation of motion ^ which satisfies the wanted initial 
conditions (|lip . More precisely, notice that 

{£2 + Xl{r) + m'a\l-\r) = and + fc^) xl{r) = . 

Hence, let us apply the operator realizing the equation of motion ([2D to the truncated series 
(Xfc — Xfe")- Using the inequalities given in the proposition, we obtain that the absolute value of 
the reminder of the latter operation is smaller then aQTn^(m^o)^"/(?^!\/2A;), which tends uniformly 
to zero over [0,to] for large n. At the same time, for every n > 0, limT-_j._oo e^'^'^Xfc (''") = 'ind 
the same holds for the first r— derivative of Xk- 

What remains to be shown is the convergence and the validity of the estimates given in the 
statement of the proposition. We proceed analyzing the n— th element of the sum. Looking at 
the recursion relations ()25p . it is a straightforward task to obtain 



\xm\ ^f/-^ rn'\xl-\t')\ dt' , \xlit)\ < j\T(t) - T{t'))a{t')m^\xt\t')\ dt' . 

The exponential bounds presented in the proposition can be obtained noticing that, W{ti, . . . , t„) 
being any totally symmetric integrable function, 

[ W{tu...,tn)dti...dtn = —.f W{tu...,tn) dti...dtn . (26) 

J0<ti<t2 - <t„<t ™' J[0,t]" 

Hence, out of (j26l) . for any V{t), 



ft ft f't\ ft 

1+ dtiV{ti)+ dtiV{ti) / dt2V{t2) + -- - =exp / dt^ 
Jo Jo Jo Jo 



V{t, 



Together with IXfcl = (2^) 2, the last observation entails both the absolute convergence of the 
series and the following bounds 

\xm\ < ( fm\T{t) - T{t'))a{t')dt'y , 

v2/c n! \Jo / 

\Xk{t)\ < ^exp /" "^aityt' , \xk{t)\ < ^exp /" m\T{t) - T{t'))a{t')dt' . 
V2k Jo k ./2k Jo 



Finally, since a{t) is growing monotonically and since r(f) —T{t') = f^, a{t") ^dt", the estimates 
presented in the proposition descend straightforwardly. 

□ 
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The previous proposition yields suitable estimates for Xklo--, 'Ai^)- We shall present the deriva- 
tion of further useful estimates also for its functional derivative Dxklo-, t, 5X](t) with respect to 
the small variation SX of X, in appendix IA.3I We shall conclude this subsection, by discussing 
some properties of the auxiliary functional $ on 23c defined as 

^X] := 27r\^{(t>^)^,^, (27) 

which corresponds to the expression (j23p multiplied by 27r^a^ and evaluated in the proper time 
t instead of the conformal time r. For that expression we have the following proposition whose 
long and technical proof can be found in appendix IA.4I 

Proposition 4.5. The following inequality holds on "Be 

mx]{t)\ <Ci{ao,c,to)t. 

Furthermore, the functional derivativ^ of ^, with respect to the small perturbation SX of X, 
satisfies the following inequality 

\D<^[X,SX]it)\ < C2{ao,c,to)to\\6X\y , 

where Ci{£to,c,tQ) and C2{£to,c,to) are uniformly bounded for tQ ^ when I and c are fixed. 
Above II • ||,^ is the norm of the Banach space SS introduced in \4-l\ 

4.5 Existence of exact solutions of the semiclassical Einstein equations out 
of initial data given at the beginning of the universe 

In this subsection, we shall present the main result of this paper, namely the existence of solutions 
of the semiclassical Einstein equations uniquely determined by some initial conditions and by the 
form of the state o^i^o in particular. Later on we shall generalize the result to other Hadamard 
states in "^(M). The proof we shall give is very similar to the Picard-Lindelof proof of existence 
of solutions of first order differential equations together with the complication that the potential 
is a functional of the solution and not simply a function. Nonetheless, we shall show that it will 
be possible to use similar methods thanks to the estimates provided by proposition 14.51 Hence, 
let us start stating the following proposition which we shall use in order to prove that there exist 
solutions of (fTB]) . 

Proposition 4.6. Fix I and let = Ito, then, if to is sufficiently small, the image of "Be under 
the map T defined in ()20p . is contained in Be- Furthermore, on 7 is a contraction. 

Proof. First of all we have to show that T maps Be to Be- In this respect, proposition 14.21 and 
the subsequent remark ensure that we can compute ^ for every element X of hence we have 
to prove that 

||T(X)-/o||,^<iy^, 
A nice introduction of these concepts can be found in [Ha82| . 
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where the norm || • ||_^ is the one introduced in definition 14.11 of B^, and /g is the center of the 
ball "Be introduced in (fT9]) . Fixing the two constants Ci = 120rrPH~'^ and C2 = (1 — c)H~'^, the 
previous inequality is equivalent to 

-C2 + (2 - c)X\t) < Ci^^m < x\t) yt G [0,to]. 

Since X'^ < t'^ ^ t^, ii is sufficiently small, this chain of inequalities is satisfied if 



\m\ < ^ 



Vt G [0, to] 



which holds true for a certain to thanks both to the bounds, satisfied by the elements of Sc, 
given in proposition 14.31 and to the estimates for found in proposition 14. 5[ 

The next step consist to show that T is a contraction on with respect to the norm of 
that is, for every Xi and X2 in 23^ 

||T(Xi) -7{X2)\\aj < C||Xi - X2\\,j , 

with < C < 1. This property descends once more from the results stated in proposition 
noticing that for every < A < 1 we can consider 

Xx := XXi + (1 - X)X2 = X2 + X[Xi- X2] , 

where Xx is in Sc. Furthermore, indicating by 6X the difference Xi — X2, we have 



< 



d_d_ 
dXdi 



7iXx)it)dX 



< sup 

Ae[o,i] 



Dj7{Xx,dX) (t) 



In order to conclude the proof we have to analyze the functional derivative of ^ in "Be- In 
this latter step, the contribution which requires some care is the one arising from the functional 
derivative of X'^a~^^, which, using both the estimates ([33]) derived in the the appendix and the 
result of proposition 14.5] can be shown to satisfy the inequality 



sup 

te[o,to] 



D 



X* 



[X,5X] (t) 



< C3to\\6x\y , 



for every X in 23c, where C3 is some constant. It is not difficult to show that similar estimates 
hold also for the other contributions to the functional derivative of Hence, there is certainly 
a sufficiently small to for which C < 1 and hence the map T acts as a contraction therein. □ 

Theorem 4.1. Fix i and let ao = ito- If to is chosen sufficiently small, it exists a unique 
solution X = H^^ of (|16p in "Be which can be constructed recursively starting from the massless 
solution Xq := Hq^ implicitly given in (|17p . In other words the sequence 

Xn := 7 (Xn-l) , 

converges in 23c to a solution of the semiclassical Einstein equations. 
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Proof. On account of proposition 14.61 t^is proof descends out of a straightforward application of 
the Banach fixed point theorem to the contraction T on Sc. □ 



Before concluding this subsection we would like to briefly comment on the result obtained 
above. First of all notice that the method to obtain the solution is constructive, in the sense 
that for every n, is closer to the solution then Furthermore, notice that if we start 

from Hq, which is a smooth function in Sc, every Hn is also smooth, because 7 maps smooth 
functions to smooth ones. Unfortunately out of the previous theorem, we cannot conclude that 
the found solution is smooth. In order to address this problem we should discuss the convergence 
of the series X„ in the Frechet spaces of smooth functions C°°[0, to], but this is out of the scope 
of the present paper. For our purpose, on account of the preceding discussion, it is enough to 
notice that we can find a smooth spacetime X„ as close as we want, in the norm of to the 
the exact solution X provided by the previous theorem. Furthermore, the regularity displayed 
by the generic element of !Bc is enough in order for equation (jl6p to be meaningful. 

As discussed at the beginning of this section we have used only part of the information present 
in the semiclassical Einstein equations, namely only the trace. What is left is completely set 
once the parameter ao is fixed in agreement with the expectation value of Tqq. Let us start 
noticing that in equation ([16]) nothing depends explicitly on oq and actually, also in {(iP')^}^^ 
as given in equation (I23p . ao does not really appear. The latter statement becomes evident, 
re-scaling everything with respect to a^, that is, using k := a^^k, r := aor, a := a^^a and 
X^(t) = ^/aQXk{T) which now satisfy the equation 



Hence, once a solution of ()16p is obtained, we can always find other solutions by choosing a 
different oq. In this sense, oq is really a free parameter that can be fixed using the information 
stored in the expectation value of Too- Actually, it can be shown that a;i^o(2oo) diverges when 
H = X~^ — )■ oo; hence, the choice the state uJi^ really provides the exact solution presented 
here, in the sense that the latter condition permits really to start the iterative construction of 
H = X^^ out of the upper brunch Hq of the massless solution ()17p . 

We conclude the present section noticing that the same result about the uniqueness and the 
existence of the solution of the semiclassical Einstein equations holds also for other choices of 
Hadamard states in '^(M) (see definition 13. 1|) provided suitable conditions are satisfied. We 
have in fact the following theorem 

Theorem 4.2. Let ua,b G '^(M) he constructed out of A and B in C^(M+), and such that, 
with respect to the auxiliary function f{k) := A{k)B{k)k'^ , the following conditions are satisfied 

a) B is rapidly decreasing, 

h) A{k)B{k) and \B{k)\'^ are contained in {W^ , kdk) , 




c) /(O) = /'(O) = 0, 
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d) f",f GL\R+,dk). 



Let Oo = ito; then, if tQ is sufficiently small, a unique solution of the semiclassical Einstein 
equations exists. 

Proof. The proof of the present theorem can be performed following the one of theorem 14. 11 For 
this reason, here, we shall discuss only the differences. Let us start noticing that what changes 
in this case is the value of (c/)'^). More precisely, instead of dealing with the map T, we have to 
deal with T' which is defined as follows 

7' {X) := 7{X) + ^^'^ 



1-H^X^ a2 

where := 27r^o^ {{(f^)u>A b ~ o) which can be rewritten as 

/■oo |- -, 

^A,B{i)= 2\B{k)\''x^Xk{t) +A{k)B{k)xJ^{t) + A{k)B{k)xkXk{t) kUk . 
Jo L ^ 

Since we already have the estimates for Xk ~ Xk given in the appendix IA.31 it is useful to divide 
^A,B in two parts, namely 

^A,B = <^A,B + ^A,B 



where ^ is constructed as ^a,b with Xk the place of Xk- Let us introduce the closed set 
formed by the elements of ^ such that 



B'^ := |x G ^ , l-td< j^X{t) < 1 + tdj . 



We have that satisfies on S^, similar inequalities as the ones stated in proposition 14.51 for 

<I>. In order to verify the last statement, we could make profitable use of the fact that B[k) is 
rapidly decreasing and that B{k)xk and A{k)xk ai^e locally square integrable with respect to 
the measure k'^dk. Hence the first bound 

\^XBit)\<Ct, VtG[0,to] 

descends straightforwardly from the inequalities in proposition 14.41 while the second one, about 
the variation of under the small perturbation 6X of X, 

mXj,{t)\<ct45x\y, vtE[o,to] 

can be proved using the results stated in the appendix IA.3i g does not satisfy similar 
inequalities, and it is for such a reason that we have introduced S^, a closed set different then 

'A,B' 



"Be- For n, there exists a constant 61 such that 



KsWi<'^i, ytG[o,to] 
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and this inequality descends straightforwardly form the properties a) and b) stated in the hy- 
potheses and from the fact that \Xk\^ constant. The variation of ^ under the small 
perturbation 6X of X involves only the variation of e'^^^'^, present in the second and third term 
From properties c) and d), it descends the following uniform decay 



< 



C 



Vt G [0,to]. 



Furthermore, operating in a similar way as in appendix IA.2t since < we obtain 

the uniform estimate 



sup 

te[o,to] 



DT{t) 



< C sup 

i'G[0,to] 



t 



\^x\ 



valid for sufficiently small io in and for some constant C . Combining these arguments we 
easily obtain 

\D^%j,(t)\<C\\8XU. VtG[0,to] 

Finally, let us notice that the results given in proposition 14.51 hold also on 23^. 

Although the behavior of ^\ ^ is worse then that of ^ or of ^a.b, the inequalities given 
above are sufficient in order to obtain the existence and uniqueness of the solution in 23^. This 
is possible because, on 23^, |A'a~^(t)| < C which is not valid for the generic element of 
Actually, for a sufficiently large d (depending on 6i in particular) and a sufficiently small to, the 
map T' introduced above, is contractive in as it can be seen adapting the proof of proposition 
14.61 to encompass the case under the present investigation. Out of this observation the proof 
of the present proposition can then be easily concluded, by applying the Banach fixed point 
theorem on 23^. □ 



5 Further considerations on the semiclassical solutions 

5.1 Expectation value of (0^) at late times 

Out of direct inspection we notice that, at late times, the massless solution (fT7|l tends to the de 
Sitter spacetime. We would like to see if such behavior persists also for massive fieldii. To this 
end notice that in |DFP08| it is argued that, in the case of a massive conformally coupled scalar 
field, if the expectation value of (j)'^ has the form f3m? + aR then, the same behavior at late times 
can be inferred. Such an estimate was obtained analyzing the adiabatic approximation for the 
state. It would be desirable to see if the assumption made in [DFPOSj for the expectation value 
of 0^ is reliable for the exact solution of the semiclassical Einstein equations. 

We start noticing that, in the case of an exact de Sitter spacetime, the expectation values of 0^ 
in the Bunch Davies ^BD78j state has exactly that behavior, and hence, there is only one choice 
of the cosmological constant which solves the semiclassical Einstein equations. Unfortunately, it 
is very difficult to compute the evolution of the expectation value of cp'^ in the exact solution we 

*See also the work |PRV08| . 
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have found before. Despite the difficulties present in the analysis of the semiclassical solution 
at finite times, we would like to stress that if a de Sitter phase exists at late times, it has to be 
exactly the one found in jPFPOSj. Such a remark originates from the fact that the expectation 
value of (fP', in an expanding flat universe {H > 0), is state independent (at least on the class of 
homogeneous pure Hadamard states) as asserted in the following theorem. 

Theorem 5.1. In an expanding H > flat FRW with smooth scale factor a which diverges in 
the limit t — )• oo, choosing a quasi-free state uja,b whose two-point function is as in (I13p and of 
Hadamard form, the expectation value of (jp' in uja,b equals the one computed on Uifi. 

Proof. The proof of this theorem can be obtained considering the form Xk- Since both uj^^b 
and o^i^o are Hadamard states, it holds for their two-point functions 



(27r)3 a{tf 



2|S(A;)r X^Xkit) + A{k)B{k) xioait) + A{k)B{k) XkXkit) dk (28) 



where x = (t,x) in cosmological coordinates and k = |k|. Because of theorem 13. 2^ B{k) is 
rapidly decreasing and A{k) of polynomial growth, hence we can pass the coinciding point limit 
under the sign of /c- integration. We aim to prove that the preceding expression vanishes in the 
limit a — >• oo. To this end we shall show that the previous k integration is bounded in time. To 
prove the last statement we have to analyze the behavior of Xk for large times. To this avail, let 
us consider the functions Xk in cosmological time {xk{t) := {xk ° ''")(0)) ^-^d let us notice that 
they have to satisfy the following equation 

Xk + Hxk + + m^^ Xk = 

where Xk and Xk indicate respectively the first and second derivative of Xk with respect to the 
cosmological time t. Hence the previous equation descends straightforwardly out of a change of 
variables in Let us introduce the positive quantity 



Q{t, k) := XkXk{t) + ( ^ + ) XkXk{t) ■ 



We would like to show that, at fixed k, Q{t,k) decreases in time. To this end notice that, since 
H>0 

Q{t, k) = -2H\xk{t)\^ - 2^H{t)\xk{t)? < 0, 



where we have indicated by Q{t, k) the §iQ{t, k). Hence, for every t > to, we have that 



Q{t,k) = Qito,k)+ [ Q{t',k)dt' < Q{to,k) . 
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Out of the form of Q{t, k) we obtain that, for every t > to, 

Q{to,k) 



\xk{tr < 



It is now possible to use the preceding inequahty to estimate the k-integral present in ()28p . and 
since both Q{tQ^k)\B{k)\'^ and Q{tQ,k)\A{k)B{k)\ are integrable in fc, we have that it exists a 
positive constant C(to), such that, for every x in the future of the surface Stg, 

\u:A,B{'P\^))-^M<t^\^))\<^ ■ 

Out of the preceding remark, we can straightforwardly evaluate the limit t — )• cxd and, since a{t) 
diverges per hypothesis therein, the difference vanishes. □ 



5.2 Variance of (T) when smeared on large volume regions 

As discussed in the introduction, the semiclassical Einstein equations are valid when the quantum 
fluctuations of the expectation value of the stress tensor are negligible. 

In this paper, in order to obtain the dynamics of the spacetime, due to the large symmetry 
present in the model under investigation, we have used only the trace of the stress tensor. Notice 
that the anomalous terms in the trace are c— numbers; hence, they have vanishing fluctuations, 
and this already means that, in the case of massless fields, the found semiclassical equation and 
its solution are meaningful at every times (not only when the curvature is small). 

Furthermore, in the case of a massive field conformally coupled with the curvature, the 
source for the fluctuations in the trace of the stress tensor can only be the expectation value of 
Unfortunately we have to notice that point-like fields like (^^(x), although represented by 
smooth functions, could have divergent variance. This is indeed the case even if the Hadamard 
regularization is employed, as can be seen noticing that in 

4>\x) (t>\y) = <P'{x)<i?{y) + 4J{(x, y)<P{x)<P{y) + 2?C(x, y)-K{x, y) (29) 

the product 'K{x, y)'K{x, y) diverges in the coinciding point limit. Above IK is the Hadamard 
singularity and -kji is the product in the extended algebra of observables =^(M), as discussed in 
[BF091 lBDF09j . Some further comments on such a product are given in appendix [ATTJ We shall 
use that expression in the computation of the variance of 0^(/), namely on a stat u it looks like 

AU<p\f)) := u:i<p\f) cl>\f)) - u{4>\f))u{ct?U)) ■ (30) 

Furthermore, if the state is quasi- free, indicating by W the smooth part of the two-point function, 
it becomes 

j [2W{x, y)W{x, y) + 4J{(x, y)W{x, y) + 2-K{x, y)'K{x, y)] f{x)f{y)dfi{x)dfi{y) 
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which, on the state we have chosen, is equal to 



where ojf q is the product of distribution ojifi • ujifi. Although q is a well defined distribution, 
due to the form of its wave front set, its coinciding point limit is not well defined and actually it 
diverges. Hence, the expectation value of (j)'^{x) has divergent fluctuations, but, luckily enough, 
we can overcome this difficulty by a suitable smearing procedure we are going to introduce. 

To this end, we could use once again the spatial symmetry and consider a smearing function 
/ni,n2 such that, for large ni and n2 it tends to have support on the whole Cauchy surface S,-, 
and it is normalized in such a way that its integral on M is equal to one. More precisely, we 
consider a compactly supported smooth function / G C^{M) centered in x,- := (t, 0), and such 
that 

f{xr) = l. [ fdfi{g) = l, 0</(x)<l. 
Jai 

Notice that, in writing the preceding expression, we are using the conformal coordinates on M. 
We can then generate the desired fni,n2 as 



X \ y/\g(ni{T' -r) +r,x/n2) 



/„,,„2 r ,x = — / -T)+r,— ^ 31 

and, as anticipated before, for large ni, /ni,n2 tends to have support on S,- while, for large n2, 
even if its spatial support becomes larger and larger, J^^ fni,n2dfJ'{g) = 1. 

We can then smear the fields in ([16]) and subsequently in pO]) with fni,n2 ^iid analyze the 
behavior of the smeared quantity for large ni and n2- In other words, notice that, since the 
trace of the Einstein tensor G is smooth and equal to —R, in the limit (ni,n2) — )• oo 

hm (5'^'^G^.,/;„„2) = -^W' 

(ni,n2)— >-oo 

and, thanks to the continuity properties of R, the result does not depend on the order in which 
the limits are taken. The analysis of the trace of (T) follows similarly. In fact, the anomalous 
part of the trace depends continuously on the coefficients of the metric, and its derivatives, up 
to the second order, while the expectation value of {(fP'{x)) has vanishing wave front set and, 
hence, it is continuous too. Thus, the equation —R = 87r(r)t^^ Q, together with its solutions, 
holds exactly in the same way for the smeared quantity. 

We proceed to analyze the variance of these quantities, and we notice immediately that the 
variances of the geometric entities present in ()16p vanish for large ni and n2, no matter in which 
order the limits are performed. We have to be more careful when analyzing the variance of the 
expectation values of {(jP') but, luckily enough, we have the following theorem which shows that 
the fluctuations of d)^ vanish if the limits are taken in a suitable order. 
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Theorem 5.2. Let fni,n2 some smearing functions constructed as in ()3ip out of some f 
with the properties stated above. Fix ni, then the variance o/ 0^(/„^^„2)? computed according to 
formula (jSOp and evaluated on the state uji^q, vanishes for large n2 in the weak sense, that is 

lim A^, o((?:)2(/„,,„2)) = . 

n2— ^oo 

Proof. First of all notice that, uJifi is of Hadamard form and hence its wave front set satisfies 
the microlocal spectrum condition. This implies that the pointwise product of oJi^ with itself, 
that we shall indicate by q, is a well defined distribution because the Hormander criterion for 
multiplication of distributions is fulfilled (the sum of the wave front set of the distribution oji^ 
with itself does not contain the zero section). Furthermore, theorem 8.2.10 of [H689] yields that 
the wave front set of the product q has to be contained in {WF{LOifi)U{0}) + {WF{ujifi)U{0}), 
where the sum is meant in the cotangent space. This implies, that the singular support of ujI q 
can only contain points connected by a lightlike geodesic. Hence, we can divide in two parts the 
distribution 



1,0 ifni,n2 ® /ni,n2) •— ^^1,0^ (/ni,n2 ® fn\,n2) + '^1,0''^ (/ni,n2 ® fni 



n2/ 



Above we have split the distribution q multiplying the test functions by a suitable partition 
of unit u + n' = l on MxM and using the linearity of oj\q. We shall now discuss how we have 
to construct u in order to have that the intersection of the singular support of wf qu' with the 
support of fni,n2® fnx,n2 is empty for every n2 and for fixed ni. Consider a compactly supported 
smooth function p on M'^ equal to one on a neighborhood of the origin. Then we shall indicate 
by u the smooth function on M x M, 

u : ((r,x); (r',x')) ^ p(x - x'), 

where we have used standard coordinates introduced above. Let us start analyzing cof qU, and 
suppose to test it on a test function h G C^{M x M) of the following form 

^((n,x); (ri,y)) := /j(ri, T2)/+(x - y)/_(x + y) , 



where ft G C^{I x /) and /+,/- € C^(M ). Furthermore, these functions are chosen in such 
a way that h is positive and, for a fixed ni, h > /ni,i (X" /ni,i- Hence, due to the translation 
invariance satisfied by cof q and the multiplication by the cutoff u, the following continuity 
condition holds true 

|^?,0«(/i)|<c[^||Z)"/t||oo)||/+||l 
\|a|<9 / 

for suitable values of C and q. Notice that, because of the localization introduced by u, C does 
not depend on the support of /_. Out of this observation we have that, for /i„ S C^{M x M) 
defined as /i„((ri, x); (n, y)) := /i((ri, x/n); (n, y/n))/n'^, 

lim ojf qU (hn) = 0. 
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The latter limit implies also that lim„2-^oo Qi^{fnx,n2®fnx,n2) — because, for fixed rii, thanks 
to the special choice of ft-, /+ and /_ and to the positivity of u^iqU, 

We have now to analyze the second contribution to i^f Q' namely lo^qu', and we notice that, 
per construction, the singular support of qu' has vanishing intersection with the support of 
fni,n2 ^ fni,n2- Hence, when tested on fni,n2 ^ /ni,n2) ^10^' has a smooth integral kernel, and, 
for every n2 and for sufficiently large ni, 

'^l.O^ {fni,n2 fni,n2) 

is represented by an ordinary integral. We can change the variable of spatial integration in such 
a way that the preceding operation can be rewritten as 

lim / {ujjQu'){Ti,T2,n2{-x-y))nlf{niT,x.)f{niT',y)dfx{Ti,T2)dxdy. 

n2-^oo J ' 

Since / has compact support and (a;f qu') is bounded, we can pass the limit under the sign of inte- 
gration. We obtain that the limit for large n2 vanishes provided that lim„2^oo(wf gu')(ri, r2, n2(x- 

y)) = o. 

Due to the spatial translational and rotational invariance, it is enough to check this last 
requirement for uji^ (without the square), for a single direction and when one of the points is in 
the spatial origin of the employed coordinate system. Hence, using the conformal coordinates, 
let us consider x := (ti,0) and y := (T2,r, 0, 0), with r >> 0. We can rewrite the two-point 
function cjifi as 

sin(A;r) 



wi,o(2;,y) = ^2 f ^ / N / Xk{Ti)xk{r2) k dk , 

where, we have made use of the fact that r is chosen to be strictly positive and the fact that (x, y) 
is not contained in the singular support of wi^o- We can now split the /c— integral on (0, oo) in 
two parts, as similarly done in the proof of proposition 14.51 namely on the interval (0, 1) and on 
(l,cx)). Since Xk{Ti)Xk{T2) is locally integrable in the measure k dk, due to the estimates given 
in proposition 14.41 and, making use of the Riemann Lebesgue theorem, we can conclude that the 
contribution on (0, 1) vanishes when r diverges. Let us proceed to analyze the A;— integral on 
the remaining interval (l,oo). Hence, we can now make use of the perturbative construction in 
for the Xfc given in (I24p and analyze the different contributions separately. First of all we 
notice that the sum of the contributions to Xk{'Ti)Xk{T2) of order equal or larger to lie in 
L^{R,kdk), and, once more, due to results stated in proposition 14.41 by the Riemann Lebesgue 
lemma, as r diverges this contribution tends to zero. We are left with the analysis of the order 
1 and m^. The former consists of a well known distribution that vanishes for large r, whereas 
the latter (order m?) can be estimated, using an argument similar to the one used to treat (j42p . 
Therefore, for large r and for some constant C, this contribution can be shown to be bounded 
by Cr~^(l -|- log(r)) which tends to zero when r diverges. □ 
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This discussion provides a new and stronger physical interpretation of the semiclassical Ein- 
stein equations on spacetime with large spatial symmetry as the cosmological one. 

6 Summary of the results and final comments 

In this paper we have rigorously constructed some solutions of the semiclassical Einstein equa- 
tions in cosmological spacetimes. All these solutions display a phase of power law inflation after 
the initial singularity. We have also seen that, if a stable de Sitter phase is asymptotically present 
in the future, then the acceleration of the latter does not depend on the particular homogenous 
pure state in which the quantum theory is evaluated, and hence it can be considered as a uni- 
versal character. Finally we have discussed the interpretation of these equations, showing that, 
when smeared on large spatial regions the equation —R = 87r(r) continues to be valid, and, in 
this case, since the fluctuations of (T) tend to vanish, it acquires an even stronger interpretation. 

We have been able to derive those results because we have selected a class of spacetimes and 
we have used an universal quantization scheme for the matter fields on that class, employing 
ideas typical of local covariance |BFV03] . Furthermore, here we have also been able to give 
unambiguously a quantum state on every element of the class of considered spacetimes, using 
and even generalizing, the results in |DMP09al IDMPOQbj . Furthermore, the employed methods 
permit to characterize the quantum states in an intrinsic way. In other words, using a lightlike 
initial surface and giving initial data thereon, we provided a prescription which can be used to 
construct states with a good ultraviolet behavior. 

There are of course many questions that are still open. In principle, similar ideas could be 
used to treat also more complicated problems like, for example, the backreaction of collapsing 
matter forming a black hole. Another interesting aspect, that should be addressed in the future, 
is the relation between the results obtained here on the vanishing of the fluctuations by smearing 
the fields on infinite volume regions and the approach typical of the so called stochastic gravity 
|HV08j . 
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A Appendix 

A.l Deformation quantization and the extended algebra of fields 

The standard procedure for quantizing the system will be to choose the symplectic space of 
the theory (§(M),cj) and then to construct the C*-algebra generated by the Weyl operator, see 
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for example |Di80j . Here we shall follow another path, that is we introduce the quantization as 
it is custom in the so called deformation quantization |BF091 IBDF09| . Let us define the off shell 
Borchers-Uhlmann algebra as 

oo 

s^{M) := 0Co°°(M)®s'^ 

n 

where only sequences with a finite number of elements are taken into account, ®s stands for 
the symmetric tensor product and the first space (n = 0) in the previous direct sum is C. The 
elements of (M) can be alternatively seen as the functionals on smooth field configurations 
if. To wit, to every F := {F^} S i/(M), we can associate 

P'if) ■= ^—AFn{xi,...,Xn),fixi) f{Xn)) , 

n 

where (p is in C°°(M) and (•, •) is the standard pairing. Out of this transformation the set (M) 
is interpreted as the set of functionals with smooth functional derivatives F^") and with only a 
finite number of non vanishing F^"'\ If we indicate by ^^"^(0) the functional derivatives of F 
we notice that F(")(0) := F„. 

The set i/(M), equipped with the pointwise product becomes a topological *— algebra with 
the complex conjugation taken as the *— operation and whose topology is induced by the one of 
the compactly supported smooth functions. Such algebra is considered to be "off shell" because 
no information about the equations of motion has been employed in its construction. At the 
end, the equation of motion needs to be implemented singling out the ideal y generated by 
elements P^^fixi, . . . , where P^,-- is the operator realizing the equation of motion applied on 
the variable Xi , we shall call this last case "on shell" . 

The quantization of {M) is realized introducing a particular star product on (M), namely 
let F and G be two elements of ^ (M) 

F*G:=^^(f("),^®'^G(")) (32) 

n 

where E is the causal propagator (the advanced minus the retarded fundamental solution). 
Furthermore, also (=e/(M),*) is a topological *-algebra with the complex conjugation as *- 
operation. 

In order to analyze the backreaction we need to deal with the stress tensor and similar 
fields. Hence, we need to introduce pointwise products of fields (like the one involved in the 
stress tensor) in the algebra of observables. Unfortunately this is usually not possible in £/{M), 
because local fields, like {f{x)5{x, y), ip{x)ip{y)) , cannot be multiplied with respect to the product 
rule given in ([32|) . Yet, since this kind of fields are needed in our discussion, we shall proceed as 
follows, first of all we isomorphically deform the algebra and then we enlarge the deformed one 
by including these and other more involved fields. The first step consists in the deformation of 
the ★ product introduced in ([32]) to the so-called Product. The latter is realized substituting 
in ()32p E (the causal propagator) with —2i'K. The new algebra (j2/(M),*j{) turns out to be 
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isomorphic to the old one; further details about this isomorphism can be found in [BF09j . After 
the deformation, we can safely enlarge the algebra (i/ (M), to {J^{M),-k'x) in such a way that 
more general local fields like the components of the stress tensor are included in (^(M), To 
be more precise, ^(M) is defined as the set of functionals over smooth field configurations, i.e., 
they are infinitely often differentiable and every F G ^(M) has symmetric functional derivatives 
F^^\0) € £'(M") (the set of distributions with compact support). Their wave front set satisfies 

(0)) n (F+ U F") = , 

where are the subsets of the cotangent space T*(M") formed by elements whose fiber 
components are future, respectively past, pointing causal or null vectors. 

We proceed now to discuss the form of the states and their interplay with the previously 
discussed deformation. Following the discussion presented in |BDF091 IBF09j . when we consider 
the enlarged algebra, we have to restrict the class of admissible states to the one satisfying 
the Hadamard condition. Furthermore, the deformation must not be visible on the expectation 
values of the original algebra. Thus, in order to preserve the expectation values of the elements 
of ^(M) under the deformation, the states on (i2/(M),*5<) are constructed out of those on 
(j# (M), by means of the push- forward induced by the isomorphism between the two algebras. 
Let us give a closer look to the effects of this push-forward on the n— point functions for the 
deformed algebra. To this end, let us fix a state u on {^Z (M), then we shall indicate by : w„ : 
the n— point function used in evaluating the elements of {£/ (M) , •k'x) and we shall call them 
regularized n— point functions. Notice that all : ujn ■ are obtained from the non regularized 
ones by the generating formula [BFOOl IHWOl) 

: Un ■■ {xi, ...,Xn) ■■= ^ , , — ( exp (l^if ® /) - Mf)] ) 

l''6f{xi)...5f{Xn) \ \2 ) ) f^Q 

where eventually the u {^p{xi) ip{xn)), needs to be substituted by ujn- Notice that, if u> 

satisfies the microlocal spectrum condition, then the all the : LVn '■ are smooth functions, hence 
the state u can be safely extended to {,^{M),-k^). Furthermore, the regularization employed in 
this paper coincides with the push forward of the state under the deformation described above, 
in fact it holds that : UJ2 := co — ^i. Finally, we stress once more that, the stress tensor T^^ and 
0^ constructed in the main text of this paper are elements of ^(M). 



A. 2 Functional derivatives of a and r 

In this subsection we shall derive some estimates, for the functionals a and r, needed to establish 
the proofs of the main theorems. Let us start by recalling the form of the scaling factor a as a 
functional of X. Fixing the initial condition at the cosmological time to and its first functional 
derivative we can write 

a[X] (t) := aoexpj-^ , Da[X,6X] (t) := a [X] (t) X-'^{t')6X{t')dt' 



35 



and, per direct inspection, we have the following estimate valid for every X in and t G [0, to]) 



\Da[X,5X] {t)\ < 



a[X] (t) 



1 1 



ao 1 



sup X-\t')t"\ ||(5X||oo < ^^\\SXU. 

i'G[0,to] *0 C 



(33) 



Furthermore, out of both the preceding inequalities and the definition of functional derivative 
as a directional derivative, it holds that 



\aiit) - a2{t)\ 



Da [Xx,dX] {t)dX 



- -r^W^^ Woo, 

to 



VtG [0,to] 



where SX = X2 — Xi and Xx = Xi + XSX and the last inequality descends from (|33|1 together 
with the proposition 14.31 

Let us continue analyzing the form of the conformal time as a functional of a and the first 
functional derivative of r o a[X] 



T[a]{t) :=ro- / -dt' , 
Jt a 

and it satisfies the following inequality 



D{Toa) [X, 5X] = Dt [a. Da [X, 5X]] 



\D{Toa)[X,5X]{t)\<2 r -{--^dt' sup \X-^{t')t'^\\\5X\\oo ■ 

Jt « V* ^0/ t'e[o,to] 

Hence, multiplying by a sufficiently large power of tm. we obtain that the inequalities. 



2tR 



t^'^D{Toa)[X,5X]{t) <C{ao,to,c) sup \X-^{t')t''^\ \\SX\\oo < ^ ||5^||oo , (34) 

t'G[0,to] «0 C 

hold in the ball "Br constructed out of c as introduced in the definition 14. 1[ 



A. 3 Bounds satisfied by Xk and its functional derivatives 

In the following we shall derive some useful estimates satisfied by the functions Xk constructed 
perturbatively in (j24p . this appendix is meant as a completion of the estimates presented in 
proposition 14.41 where the convergence of (j24p was proven. All the forthcoming formulas are 
meant to be valid for every spacetime that arises from an element X of 23^ introduced in definition 
14. H along the lines of proposition l4.1l (Notice that they holds also on "B'^ introduced in the proof 
of theorem 14. 2p . For the same reason t G [0, to] everywhere. For our later purpose, in order to 
shorten some formulas, it is useful to introduce the following two functionals 



FMit) := l^in^dt' , F2[a]{t) := 



(r(t) - T{t'))a{t')m'^dt' 



(35) 
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that satisfy |Fi[a](t)| < -^t'^rn? and |F2[a](t)| < m?t'^. Following a procedure similar to the 
one presented in the proof of proposition 14.41 can derive the following estimates 



\xnm\ < ^(Fi/2[a](t))"+^expFi/2[a2](t) . (36) 



1 

2k 

In the remaining part of this appendix we shall analyze the functional derivatives of Xk ~ Xk 
with respect to X G "Be, more precisely we shall compute 

D{xk - xl) [X, 5X] = Daixk - xl) [a, r, Da [X, 5X]] + D^ixk - xl) [a, r, Dt [X, 5X]] . 

Let us start from Da{xk — t, Da[X, (5X]], which, exploiting the properties of the recursive 

series ([M]) . the form of Da[X,5X] and the identity ([26]) . turns out to satisfy the following 
inequality 

\Da{xk-xl)[a,T,Da{X,5X)]{t)\ <—j=mH\\Da\U 

where || • ||oo is the uniform norm on C[0,to]) while Daxi ~ ^- '^^^ second contribution is more 
laborious; let us consider the recursive sum and let us concentrate on the contribution to the 
n-th order which looks like 

f e-M^") fr sin(fc(r(t,-i)-r(tO)) , 

/ — r^li ; m a{ti)dti . . . dtn (37) 

Jo<t,,<t„-i-<t y2k .^-^ K. 

where, only in the previous formula, to corresponds to the time t. The action of the r-functional 
derivatives on (I37p is twofold: On the one hand, we have to derive every factor of the form 
sin(2A:(r — r')) and, on the other hand, we have to derive r in x^- Let us first of all consider 
the functional derivatives of sin(2/c(r — r')) and, using (f26|l . we obtain that the n— th order is 
bounded by 

, \„ I- F!'z\t) ■ [ rr?— ( ^] ' a{t')dt' ■ \\6X\\oo 
(n-1)! V2fc ^ Wo cao \t' J ^ ^ " "~ 

where || • ||oo is the uniform norm on the elements of ^ and where we have used the following 
estimate descending from ([M]) 

\Dt[X,5X] ih)-DT[X,5X] (t2)| < — f^V \\6X\\oo, h > t2 . 

The absolute value of the functional derivative of e^^'^ in ()37p is also bounded by a similar 
expression 



F^-\t). [ m'^(^-^Y a{t')dt'.\\6Xl 
' Jo c ao \i / 



(n-l)!V2fc 

Both these last two estimates can be summed in order to give an exponential. Hence 

, „ , exp Fi /oft) 6 m? 

DAxk - xl) [x,dx] (t) < ' I^^' t. — -toWmU 

V2A; 2c - 1 
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and, together with the a-functional derivative, we obtain 



\D{xk-xl)[X,5X] {t)\ < 



1 ao 1 6 m 



2 1 



+ 



expFi/2(t) 



2k 



m\\\5X\ 



k to c2 ' 2c- 1_ 

A further inequahty we would hke to present in this subsection is the fohowing 
{Xk-xl)[X] {t).D'7AX,SX] (t) 



< 



expFi/2(i) 2 m" 



2k 



2c -I 



to\\6X\ 



(38) 



(39) 



which can be shown to hold in a similar way as the bound found above for the functional 
derivative of e^^'^ . We conclude this subsection with two inequalities, valid for n > 0, that can 
be shown to hold in a similar way as before: 



I^xr [X,SX] it)\<F^,it) 



expFi/2(t) 



2k 



2 1 



1 ao 1 6 m 
kli^'^ ^ 2c -1 



and 



X>''[X]{t).Dxl[X,5X]{t) <F^,,{t) 



expFi/2(t) 2 



m 



2k 



2c -I 



m^toll^Xl 



tQ\\5X\ 



(40) 



(41) 



A. 4 Proof of Proposition 14.51 

Proof. The proof is done dividing the expression for <I> in a finite number of parts and showing 
that every contribution separately satisfies the desired inequalities with respect to some constants 
that have the properties stated in proposition 14. 5[ In the next, since most of the computations 
are rather long but straightforward, we shall concentrate only on the key points, and we shall 
make extensive use of all the estimates derived in the previous appendix. Throughout this proof 
we shall make use of proposition 14.11 in order to associate a spacetime to an element X of Sc. 
According to the same proposition, we shall indicate as H, and the t which appears in 
some estimates has to be thought as being contained in the interval [0, io]- Let us start dividing 
the integral over k in two parts 



^\X] 



dkk^ 



XkXk 



xlxl 



+ 



dkv 



XkXk 



xlxl + 



4fc3 



m a log(am) , 



where the regularization is needed only for large A:, whereas X) o ^ind r are functionals of X 
as discussed previously. Notice that, the zeroth order in m^, present in {(\?') through X^Xfc) 
vanish due both to the discussion about the regularization freedom and to the requirement that 
Minkowksi spacetime, with respect to the Minkwoksi vacuum, is a solution of the problem. For 
this reason we have subtracted it in (j27p and we shall not consider it anymore. Let us start 
considering the first part of the integral, to which we shall refer to as the infrared part and we 
shall indicate it by <I>^. In this respect we shall write 



XkXk XfcXfc 



Xk - Xfc ) [Xk - Xfc) + [xk - Xfc ) Xfe + Xfc [Xk - Xfc) 
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and, hence, using the estimates ([36]) for Xk and the one for the functional derivatives ([38]) and 
([39]) . with F2 as in ([35]) . we obtain 



\^\t)\ < J [F2(t)2exp(2-F2(t)) + 2F2(t)expF2(t)] <\ [mH'^ + 2mh'^) e 



2x2 



where F2{t) < rn^t^, as it can be seen by direct inspection of the definition of both F2{t) and r, 
and 



D^^ [X, 5X] {t)\<\ [F2{t) expF2(t) + 1] 



ao 1 3 m? "I 

■0 



in c2 ' 2c - 1 



+ 



1 m"' 
2c- 1 



expF2(t) m^toll^^Xll, 



This conclude the analysis of the contribution of the lower energies to the value of we can 
now proceed to analyze the higher frequencies which we shall indicate as . We shall split 
it in powers of rm? and discuss the first three powers separately. Afterwards, before summing 
everything together, we shall analyze the remainder. 

I) order. Let us start considering the first order in and its contribution to <I>^ in particular. 



$f[a,T](0 : = 



4A;3 



dk — m^a^(t) log(a(t)m) 



which can be expanded as 



/■oo 2 ft Q 2 

$f[a,r](t) := / dk — cos {2k{T - T'))-^{t') dt' - m^a{tf log {m a{t)) 
Ji 4A; Jq dt 

where both a and r have to be thought as functionals of X and where we have written the 
perturbative series (I24p using t, instead of r, as the time variable. Furthermore, r has to be 
intended as r[X](t) while r' as r[X](i'). Let us rewrite it in the following form 

dk^ smi2k{T-T'))-l^a-^j (0 dt'-m\{tf log {ma{t)) 

and, using the fact that the t derivative of X is positive, we can obtain the following bound 

<m'^a^{t)X~^{t) + m'^a{t)t + m'^a'^{t)log{ma{t)) . 
Let us proceed to analyze the functional derivative 

L>$f [a, T, 6X] = Da<^^ [a, r, Da[X, 6X]] + L>^$^ [a, r, Dt[X, 6X]], 

and let us analyze both term separately, first of all, setting F := Da[X, 6X], 

t 

sin (2/c(t - r')) F{8a^H^ + Aa'^H) + SFa^H + 2a'^F {t') dt' dk- 
-2m^a{t)Flog {ma{t)) - m^a{t)F , 



8A;2 
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where as usual H = and the dots indicates the derivatives with respect to the cosmological 
time t. Using the estimates present in proposition 14.31 and ([33|) . it is now an easy task to obtain 
the desired bounds. Notice in particular that the second order time derivative of F multiplied by 
can be estimated by (the norm of as it can be verified directly from the definition 

of a[X]. Let us proceed with the second functional derivative, and now setting G := Dt[X,6X] 
we obtain 

Dr^^[a, T, G]{t) := ^ / cos (2A;(t - r')) [G{t) - G{t')] ka^H^ + 2a^H] {t') dt' dk . 
Ji 4fc Jq L J 

We could in principle integrate by parts in order to obtain an extra l/{2k) factor. Unfortunately 
this would not help much, since, due to the extra i-derivative the left hand side would depend 
on H and we do not have control on the second f-derivative of H in "Be- At the same time we 
have to remember that is integrable in k on [1, oo) for every x > and that the fc-integral 

present above needs to be understood in the distributional sense, hence taken with the following 
e-prescription 

p OO 2 pt 

Dr<^^[a,T,G](t') := lim / e~''''^ / cos (2 A: (r - r')) \G(t) - G(t')] (6a^ + 2a^ H)(t') dt' dk . 

e=0+ J I 4k Jq 

Since [G{t) - G{t')] {6a^H^ + 2a^H) is ^^(0,*), at fixed e, the absolute value of the integrand 
in dk A dt is integrable. We can thus switch the order of integration to obtain 

Dr<^V[a,T,G](t) ■.= m'^ \G(t) - G(t')] (6a^H^ + 2a^H)(t') lim / -— cos (2k(T - t')) dk dt' . 

Jo ^=o+Ji 4k 

(42) 

We are now ready to perform the integral in k and to notice that the result is a smooth function 
everywhere in r — r' but in where a logarithmic divergence in r — r' appears when the limit 
e = 0"*" is computed. Let us call L(r — r') the result of the ^-integration in the limit e = 0. 
Notice that L{T{t) — t' {t')) / a{t') is absolutely integrable for t' in \t — 6, t] and, moreover, we can 
adjust 6 in such a way both that this integral is equal to 1 and that |L(r(t) — T'{t'))\ < 1 for t' 
in {0,t- 5). Hence, since a(t') [G{t) - G{t')] {6a^H^ + 2a?H){t') is uniformly bounded for t' in 
[0,t], we can split the t'-integral in (|42p in two parts, namely over [0,t — S\ and over \t — 6,t\. 
Both can be uniformly estimated and finally the following absolute bound for IDt-^j^] can be 
established 



\Dr^^{t)\ < sup \a{t') [G{t) - G{t')] I \\6a^H^+2a^H\\^+ [ [G{t) - G{t')] {6a^H^ + 2a^H) 
t'e[o,i] Jo 

Combining all the estimates, we obtain the desired behavior, namely that the absolute value of 
the functional derivatives is controlled by Cto||5^||- 

II) order. The computation at the second order is more involved due to a pair of integrals which 
needs to be addressed, but at the same time there are no subtleties with the extra derivatives 



dt' 
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which could necessitate control on the derivatives of the variation higher than the first one. We 
start reminding the form of this term 



xlXk [a, T-] + xlxi [«' + XkXi [a, t 



A -.A 



dk , 



which can be more explicitly rewritten as 



$^[a,r](t) := / dk — 



dk 



m 



4A:2 



/ sin(2A:(r-r'))^(t') / a{t")dt"dt' + / dk 
Jo Ot Jt' Ji 



1 



sin{2k{T-T'))^{t')dt' 



m 



4 rt 



sm{2k{T -T'))a^{t')dt' . 

(43) 



Out of such expression, an estimate satisfied by |$2^| can be easily derived and it looks like 



Let us proceed to the analysis of the functional derivatives of ^ and let us notice that the 
analysis of Da^2 can also be easily done yielding 



\Da^^[a,T,Da]{t)\ < — [Ga^H^t'^ + Ht'^ + ta^] ||L>a||oo + — [2a'^ Ht"^ + at'^] a{Dc 

8 8 

A little bit more laborious is the handling of the r-directional derivative, because, operating as 
above, an extra factor 2k appears in the formula. Furthermore, it could make the A;-integral 
divergent when the trigonometric functions are replaced with one in the maximization procedure 
of the second, third and fourth factor in ()43p . In order to avoid this problem, we can integrate 
by parts after performing the functional derivative, along the following lines: 



m 



DrJ^ sm{2k{r-T')) — dt' = s\n{2k{r - r'))- a(—\D{T-T') 



dt'. 



The price to be paid is in the extra time derivative, we have to deal with. In fact, proceeding 
in that way we obtain the following estimate 



\Dr^^[a,T,DT\{t)\<Ci I a\D{T - T')\dt' + C2 [ a 

Jo Jo 

It is now a straightforward task to obtain the desired inequalities. 

Rest2) The contribution of the order larger than 2 in m? is 



dt' . 



^'i,[a,r]:= / dk k' xlxl + xlxl + ixl + xl) xf + xfixl + xl) + xV xV 
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Hence, employing the inequality ()36p , and the definition of -Fi in (|35p , we obtain the bounds 



'A;2 



to 



1 + ( 1 + ^^/] exp(Fi[a]) + ^%'expm[a]) 



where Fi{t) < ^j^t^, more explicitly we can rewrite it as 



|^>2[a>T"](t)| < 3 ^m^^t^^ (1 + 771^00*0) exp {2m^aoto) . 

We do not analyze the functional derivative of this contribution, because, due to the extra k 
appearing in the functional derivative of the trigonometric functions with respect to r, we have 
to further split it. 

Ill) order The contribution at third order in is 



and its explicit form is 



dk k^ 



Xlxl [a, t] + xlxl [a, r] + xUl [a, r] + xlxl [«, r] 



$3^[a,r](t) : = 



dk- 



A;2 



/ sm{k{T - Ti)) sin(A;(Ti - T2)) sin(/c(r2 - t^)) cos(/c(t - Ts))a{ti)a{t2)a{ts) dtidt2dt3+ 

J 0<tz<t2<ti<t 

+ / / sin(A;(T - ti)) sin(/c(r - T2)) sva.{k{T2 - t^)) cos(A;(ti - T^))a{ti)a{t2)a{t-i) dtidt2dt3 
Jo J0<ta<t2<t 

(44) 



An estimate of the form 



D^'i{t)\<C{t^,ao,c)\\5X\ 



holds for the first functional derivative. Notice that, in the derivation of the preceding estimate 
starting from (j44p , a problem could in principle arise when the r functional directional derivative 
is computed because in such a case an extra k factor appears and the k integral is logarithmic 
divergent. Hence, in order to obtain a finite /c-integration, before performing the functional 
derivative, we have to rewrite the product of the trigonometric functions as a sum of products 
of trigonometric functions depending only on r — ti, ti — r2, T2 — Ts for the first line in (I44p and 
on r — Ti, T — T2, T2 — T3 for the second one. Subsequently we can safely perform a couple of 
integration by parts. The outcome is that the result of the t-integrations is of order l/{2k) in k 
and this extra k permits to obtain finite bounds for the ^-integration also for the r-functional 
derivative. The price to be paid is that, afterwards, in the ti integrations the derivatives of 
a{t) might appear. This also holds for the a directional functional derivative which has to be 
carefully computed. More precisely we obtain 

ft 



dk 



^ [im'^a'^H + 6ta^] [ a\Di 
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and 



\Da'^3it)\ < 



/oo "I r roD "I 

dk^ [SAHa^t^ + 6a^t^]\\Da\\^+ J dk^ [28a'^t^] a{Da) 



out of which it is straightforward to obtain the wanted inequahty. 
Rests) The contribution of the order larger than 3 in tv? is 



dk k^ 



XkXk +Xk Xk + XkXk +Xk Xk + Xk Xk 



and, using (j^T]) and ([ICTj) together with ([55]) . we can find a constant C"(ao,to) c), with the desired 
properties, i.e.. 



\D^%{t)\<C'{ao,to,c)\\5X\ 



yte [0,to] 



After a rather long but straightforward computation, we are ready to collect all the estimates 
presented above hence obtaining that 



Mt)\<C{ao,c,to)t, \D^[X,6X]{t)\ < C{ao,c,to)\\5Xl 



Vt G [0,to] , 



where C{tQ£, c, to) remains bounded for small to- The thesis of the proposition can be obtained 
noticing that ||(5X||oo < ^oll'^-'^ll,^ in particular. □ 
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